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Abstract 

This  investigation  is  concerned  with  the  development  of  a  two-body  distribution  func¬ 
tion  in  a  plasma  for  use  in  a  kinetic  equation  for  the  one-body  distribution  function.  The 
kinetic  equation  is  obtained  for  a  uniform  plasma  for  those  circumstances  in  which  the 
time  dependence  of  higher-order  distribution  functions  can  be  assumed  to  occur  within 
a  functional  dependence  on  the  one-particle  distribution  function.  The  conditions  of 
validity  for  this  functional-dependence  assumption  are  discussed.  The  resulting  inter¬ 
action  term  is  new  in  the  sense  that  it  contains  no  divergent  integrals  requiring  cutoffs, 

auld  it  may  be  considered  accurate  to  first  order  in  (e  /kTkjj).  The  interaction  term  is 

composed  of  two  parts.  The  first  is  a  Boltzmann  collision  integral  with  a  Debye -shielded 
interaction.  The  second  term  is  due  to  the  deviation  of  the  shielding  cloud  from  a  Debye 
shield  and  is  the  Fokker-Planck  form,  the  coefficients  of  which  are  finite  and  well- 
behaved.  Because  of  its  form,  with  a  convergent  collision  integral  and  convergent 
Fokker-Planck  coefficients,  the  solution  may  be  considered  as  a  joining  of  the  previous 
solutions  to  this  problem. 
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I.  THE  NATURE  OF  THE  PROBLEM. 

1.1  INTRODUCTION 

In  the  description  of  the  evolution  of  a  fluid,  the  interplay  of  the  properties  of  the 
interaction  and  the  density  determine  whether  the  fluid  is  a  liquid,  gas  or  plasma.  If 
at  the  interparticle  distance  the  motions  of  the  particles  are  strongly  affected  by  the 
forces  that  are  due  to  other  particles,  we  consider  the  fluid  a  liquid.  Forces  that  are 
weak  at  the  interparticle  distance  will  give  rise  to  fluids  described  as  gases  or  plasmas. 
If  we  consider  forces  that  fall  off  at  large  distances  to  be  r”“,  we  may  distinguish  gases 
and  plasmas  as  a  >  2  and  a  ^  2,  respectively.  For  gases  the  effect  of  the  volume  ele¬ 
ment  in  a  summation  over  particles  is  not  enough  to  counteract  the  r  °-decay  of  the 
force.  The  summation  will  receive  its  greatest  contribution  from  near  particles.  In 
a  gas  of  this  type,  for  a  density  that  is  sufficiently  low  that  particles  interact  only 

infrequently,  we  are  led  to  the  concept  of  a  two-body  collision  and  a  Boltzmann  gas. 

2 

For  a  <  2,  the  r  -increase  of  the  volume  element  is  sufficient  to  increase  the  effect 

of  the  more  distant  particles.  In  particular,  for  Coulomb  interactions  in  a  plasma  the 

_2 

volume  element  just  offsets  the  r  -decay,  and  the  concept  of  a  collision  becomes  vague, 
since  many  particles  interact  at  once. 

Kirkwood  and  Poirier  show  that  for  a  plasma  in  equilibrium  the  effect  of  the  Coulomb 
potential  is  modified  by  the  screening  of  other  particles  and  an  effective  potential  is 

well  approximated  by  the  Debye  potential, 

for  n  particles  per  unit  volume,  becomes 

a  new  range  of  interaction.  As  we  shall  see  (sec.  1 .2),  plasmas  of  interest  will  be  of 
such  a  density  that  they  will  have  many  particles  within  a  radius  nk^  is  a  large  num¬ 
ber.  We  are  still  faced  with  many  particles  interacting  at  once. 

Because  the  information  that  is  eventually  desired  about  a  gas  will  not  depend  on  the 
detailed  motion  of  particular  particles,  but  rather  on  an  average  over  many  particles, 
it  is  convenient  to  introduce  distribution  functions .  Thus  we  introduce  the  distribution 
functions 

F^(x^;t),  F2(x^,X2;t),  ....  Fg(x^,  . , .  x^n). 

where  the  subscript  denotes  the  number  of  particles  in  the  argument;  t,  the  time;  and 
X.,  the  6-dimensional  vector  {p.,q.}.  We  shall  assume  that  the  F  are  invariant  under 
an  interchange  of  particles,  so  that  the  particular  x.  appearing  in  the  arguments  are 
simply  labels.  The  F  may  be  defined  as  follows:  The  function — -  F  (x,,  . . .  x  ;t)  dx, 

S  SI  S  1 

. . .  dx  is  the  probability  at  time  t  that  the  particles  1  . . .  s  will  be  found,  respectively, 
s 

at  X, ,  ...  X  within  dx,  . . .  dx  .  We  shall  take  the  relation  between  the  distribution 
I  s  Is 

functions  to  be 

F^(Xi,  .  . .  x^;t)  =  ^  j  F^(xj . .  .x^,  t)  dx^  . . .  dxg^^  k  >  s. 


established.  The  effective  potential  can  be 
•p  e  The  Debye  length,  kj^ 
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where  the  volume  V  that  is  available  to  the  gas  is  inserted  in  both  cases  to  allow  a 
smooth  passage  to  infinite  volume.  For  each  dx^,  the  integration  is  over  all  momentum 
space  and  the  volume  V. 

There  is  one  more  concept  that  will  be  of  interest  to  us.  If  the  motion  of  the  s  par¬ 
ticles  can  be  considered  statistically  independent,  we  have 


Fs<"l' 


(1) 


In  this  report  we  shall  suppress  the  time  variable  t  when  it  is  not  important  to  the 
relation  considered.  We  shall  refer  to  motion  for  which  the  property  (1)  holds  as  uncor¬ 
related  motion.  Motion  for  which  (1)  does  not  hold  will  be  referred  to  as  correlated 
motion. 

The  distribution  functions  have  final  interpretation  when  used  to  evaluate  the  average 
or  expectation  value  of  some  mechanical  property  A  that  depends  on  s  particles .  We 
interpret  this  property  as 


As(Xi. 


.X  ) 

s' 


F  (x,, 
s'  1 


Xg)  dXj 


dx. 


(2) 


where  we  have  used  the  fact  that  the  distribution  functions  are  normalized  as 

F^(Xj)  dXj  =  1.  Obviously,  we  are  primarily  interested  in  s  =  1  and  s  =  2. 

Upon  introducing  the  distribution  functions,  we  are  interested  in  their  equations  of 
motion  rather  than  those  of  individual  particles,  although  the  two  are  closely  related. 

One  of  the  first  successes  in  this  direction  was  Boltzmann's  equation  for  the  one-particle 
distribution  function  for  dilute  gas  of  short-range  potential. 


»F,(x  )  ^ 

8t  ^  ''l 


8F^(x^) 

8qi 


=  1  ^(p^)  Fj(p2,q2)}d?2adad4.. 


(3) 

Here,  a  is  the  collision  impact  parameter  and  Pj  and  p^  are  the  momenta  that  the  par¬ 
ticles  must  have  had  before  the  collision,  given  that  their  coordinates  are  now  Xj  and 
x^-  We  shall  consider  no  external  force. 

It  may  be  recalled  that  Boltzmann's  derivation  depended  upon  a  long  free  path  between 

relatively  quick  collisions.  As  an  approach  to  the  problem  of  plasmas,  this  equation 

Z 

with  a  modified  Coulomb  potential  was  used  by  Spitzer  and  Harm.  The  collision  integral 
_2 

with  a  straight  r  -force  diverges  at  long  distances.  Using  the  known  fact  that  at  equi¬ 
librium  the  effective  interaction  is  the  Debye  potential,  they  cut  the  integral  off  at  X^,. 

3  ^4 

Another  approach  was  originally  proposed  by  Vlasov  and  solved  in  detail  by  Landau. 
This  emphasized  the  Coulomb  nature  of  the  interactions  by  considering  the  force  on  a 


2 


The 


charge  to  be  given  by  the  gradient  of  a  potential  whose  source  is  J  F^(Xj)  dp^ . 
equation  for  becomes 


8Fj(Xj) 

5t 


8Fj(x^) 


^8F^  ^  n  F^(x^)  dx^ 

""  BV,  '  8qi  Iv^l 


(4) 


As  we  shall  see  later,  this  equation  is  a  first  step  but  does  not  include  or  account  for 
particle  correlations. 

A  third  approach  to  the  problem  may  be  made  through  a  Fokker-Planck  type  of  equa¬ 
tion  for  Fp^ 


9F^(Xj) 

at 


+  V, 


aF^(x^) 


a 

av 


^D°Fj(Xj)  + 


avj  J 


(5) 


where  D°  and  are  functionals  of  F^ .  The  form  of  this  equation  is  derived  by  assuming 
that  the  particles  undergo  a  large  number  of  small  deflections  —  a  condition  violated  by 
Coulomb  particles  undergoing  a  close  collision  with  large  momentum  transfer. 

Before  the  Fokker-Planck  equation  can  be  of  use,  we  must  obtain  the  coefficients 
from  considerations  of  the  interactions.  One  method,  described  by  Allis, ^  is  to  expand 
the  Boltzmann  collision  integral  in  terms  of  small  deflections  and  obtain  the  Fokker- 
Planck  form  of  equation.  This,  however,  necessitates  using  the  collision  integral  in  the 
region  in  which  its  accuracy  is  most  suspect  —  long-range  interactions  or  grazing  col- 
lisions.  Another  method  described  later  uses  the  two-body  distribution  function 
under  the  assumption  that  the  two  bodies  never  get  close  together  or  experience  an  inter- 

7 

action  that  is  strong  compared  with  their  kinetic  energy.  As  discussed  by  Balescu,  the 
coefficients  of  this  equation  have  great  intuitive  appeal.  However,  integrals  in  them 
diverge  at  short  distances  because  of  a  violation  of  the  initial  assumption;  in  this  report 
we  shall  resolve  this  divergence. 

One  more  fact  should  be  noticed  about  these  three  equations.  On  the  right-hand  side 
of  the  Boltzmann,  Vlasov,  and  Fokker-Planck  equations  we  have  functionals  of  F^  only. 
The  future  of  F^  is  determined  by  its  present  value,  but  not  obviously;  we  might,  for 
example,  find  that  the  equations  for  F^,  F^,  and  so  forth,  are  all  interrelated  in  a  set 
of  simultaneous  equations.  As  it  is,  the  equation  for  F^  and  the  higher-order  distribu¬ 
tion  functions  must,  in  some  sense,  be  trivial  in  the  time  variable  so  that  they  can  be 
solved  immediately  with  the  result  of  a  single  equation  for  Fj .  This  occurrence  will, 
be  discussed  in  great  detail  in  section  2.1. 

The  purpose  of  this  report  is  to  derive  a  kinetic  equation  for  F^ .  We  shall  derive 
this  equation  by  starting  with  the  general  Liouville  equation  for  the  plasma  and  then 
examining  the  circumstances  under  which  a  kinetic  equation  for  F^  can  be  assumed  to 
exist.  For  these  circumstances  we  shall  obtain  an  equation  for  F,  to  first  order  in  the 
small  parameter  (e  /kTX.j^).  This  will  be  an  improvement  upon  the  attempts  mentioned 


3 


above  in  the  sense  that  it  will  contain  no  divergences  in  those  terms  corre¬ 
sponding  to  interactions  with  other  particles,  that  is,  the  right-hand  sides  of 
(3)-(5). 


1.2  STATEMENT  OF  THE  PROBLEM 


8-10 

Several  authors  ~  have  used  the  Liouville  equation  for  the  distribution  function 
in  the  phase  space  of  6-N  dimensions  as  a  starting  point  for  the  discussion  of  the  evolu¬ 
tion  of  a  gas  of  N  particles  in  a  volume  V.  This  distribution  function,  Dj^(Xj,  . .  .  Xgit). 
is  assumed  to  be  symmetric  under  the  interchange  of  any  pair  of  particles.  The 
Liouville  equation  is  then 


30. 


N 


3t 


N 


KJ 

Here,  the  brackets  are  Poisson  brackets, 

[F,G]a  f  /iZ.iG.M.iGl 

itlls^i  »Pi  ^Pi 

and  the  kinetic  and  potential  energies  are 


(6) 


(7) 


i=l 


and 


4- 


ij 


For  our  purposes  throughout  this  report  we  shall  assume  only  Coulomb  interactions  and 
identical  electrons  of  mass  m  and  charge  e  imbedded  in  a  uniform  background 
of  opposite  charge.  The  net  charge  will  be  taken  to  be  neutral.  The  points  to 
be  investigated  may  be  studied  with  this  idealized  model,  without  the  compli¬ 
cation  of  different  particles. 

We  define  the  reduced  distribution  function  for  s  particles 


Fs<"l' 


3'^^  ~  „N-s  J  ‘^s+l 


dx 


N' 


(8) 


N-s 

where  the  factor  V  is  introduced  to  allow  a  transition  to  infinite  volume. 
We  shall  also  have  use  for  the  following  identities: 


4 


V 


"aT'^s+i 


N 


dx 


s+1 


<^N  •  [«=•' 


r  ^ 


LKJ 

'  S 

i<J 


(9) 


dx 


s+1 


.  dx 


N 


,N-s 


N  s 


dx 


N' 


j=s+l  i=l 


By  using  the  symmetry  of  Dj^  under  interchange  of  particles,  the  sum  over  j  consists 
of  identical  terms  so  that  the  last  identity  may  be  written 


N 


}<i 


N 


dx 


s+1 


.  dx 


N 


”  s 

rs 

1 

)  4/..;F 

^  '^is+l’^s+l 

dx  . ,  , 
s+1 

,i<J 

Li=i 

Integration  of  (6)  over 

dx^+1  ... 

dXj^  and  use  of  (9)  produce 

8F 

r  “1 

1 

■fcT  O 

*  s 

- 

s 

8t  ■ 

+ 

L  s  sj 

Z 

I 

“^is+l  '^s+1 

_i<j 

i=l 

(10) 


dx 


s+1 ' 


(11) 


Throughout  this  report  the  argument  of  the  function  is  suppressed  when  the  meaning 
is  clear. 

For  fixed  s  we  may  pass  to  the  limit  of  infinite  volume  and  infinite  number  so  that 

N  N  —  s 

remains  constant.  This  replaces  — with  n.  For  spacially  nonuniform  plasmas, 

''  ''  N 

n  is  not  a  density  in  quite  the  usual  sense,  but  is  the  limit  y-. 

The  introduction  of  n  gives  (11)  its  final  form.  The  meanings  of  its  terms  are  clear. 

The  left-hand  side  and  the  first  two  terms  on  the  right-hand  side  constitute  the  Liouville 

equation  including  interactions  for  the  s  particles  under  consideration.  The  integral 

term  represents  the  contribution  to  the  rate  of  change  of  F  which  is  due  to  the  inter- 

s 

actions  with  the  rest  of  the  particles. 

In  order  to  approach  the  analysis  of  the  equations  (11),  we  shall  estimate  the  size  of 
the  terms .  It  will  be  found  that  under  certain  circumstances  one  term  is  small  and  thus 
gives  rise  to  the  possibility  of  a  perturbation  expansion.  We  assume  a  plasma  that  is 
near  enough  to  equilibrium  that  we  may  define  a  shielding  distance  =‘^ kT/4Trne^,  a 
plasma  frequency  w  =^4irne^/m,  and  a  characteristic  velocity  v  =  By  using 
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these  as  units,  the  sizes  of  the  various  integrals  and  derivatives  can  be  estimated.  The 
change  of  coordinates  t  =  cj  t,  p  =  r/X.^,  and  u  =  v/v  in  (11)  will  make,  for  the  average 
particle,  each  of  the  integrals  and  derivatives  of  order  one.  The  size  of  each  term  will 
be  given  by  the  following  coefficients 

2  2  2 

Wp,  e  /mvb  ,  ne  X.jj/mv 

which  are  the  ratios 


1,  1,  l/nXjjb^,  1. 


Here,  b  is  the  distance  between  the  two  particles  under  consideration.  The  last  three 
terms  are  the  three  sums  in  (11),  and  the  coefficients  above  simply  represent  their 
magnitudes.  In  particular,  in  the  third  term  there  will  be  a  different  b  for  the  sep¬ 
aration  of  each  pair  of  particles;  the  term  for  each  pair  is  subject  to  the  analysis 
given  below. 

A  requirement  that  the  third  term  be  small  implies  that  b  »\/e^\y.,/kT.  The 

8-10  ^ 

solution  of  (11)  has  been  discussed  by  several  authors  '  for  the  situations  that 
satisfy  this  requirement  for  all  pairs.  The  process  has  been  to  assume  that  all 
of  the  s -particle  interaction  terms  are  small  and  to  assign  to  them  an  expansion 
parameter  g  that  is  later  set  to  one.  The  terms  designated  by  g  will  be  of  order 
(e  /kT\jj)  .  Then  is  expanded  in  g  as  a  perturbation  expansion.  This  results  in 
zero  and  first  order  in  g. 


flF- 

IT 


[«s'Fs]  "  I 


'*'is+l'^s+l 


i=l 


dx 


s+1 


and 


aF 

It 


1 

s 


’  s 

”  s 

+ 

y  i|>..;F° 

y  '^is+l'^s+l 

i=l 

dx 


s+1  ■ 


(12) 


(13) 


The  derivation  and  solution  of  these  equations  will  be  discussed;  of  more  interest 

to  us  now  is  the  general  form  of  the  solution.  We  shall  see  that  for  a  uniform  plasma 

F°  and  F^  may  be  reduced  to 
s  s 


F°  =  fr  Fi(Pi) 

®  i=i  ^  ^ 

(14) 

s 

i<j  ’  ^ 

Here,  F^  is  the  solution  of  the  equation 


6 


8F2(x,,x,) 


-  jH2;F2(Xj,  x^)  + 

^  ■*■  ['^23'^1^P2^  ^2^*1’ *3^]}’  ‘^3’ 


(15) 


and  to  first  order  in  g 


9Fl(Pl) 

8t 


(16) 


Equation  15  has  been  solved  under  the  adiabatic  hypothesis  defined  in  section  2.1.  When 
this  solution  is  substituted  in  (I6),  the  final  equation  becomes 


1^3  4 

loir  ne 


,dk  k 


6(k-(Vj-V2)) 


m(2ir)' 


8V, 


k2  + 


2  C 

(A)  \ 

P  J+ie 


aF,(P2) 

^  ■  8p, 


dp. 


k  •  (Vj-V^) 


4-^1  Fi(Pi)  F,(P2). 

8Vi  dVj  ^  ‘  ‘ 


where,  under  the  adiabatic  hypothesis. 


(17) 


(18) 


For  definiteness,  we  reserve  a  special  symbol  for  this  function  and  refer  to  it  as  the 
large-separation  solution. 

In  (17)  k  is  the  Fourier -transform  variable  corresponding  to  (rj-r^) .  This  inte¬ 
gration  diverges  logarithmically  for  large  values  of  |k|  corresponding  to  small  values 
of  I  Tj-r^  (.  As  mentioned  previously,  this  divergence  occurs  because  we  assumed  that 
v|(j2  Is  small  in  the  derivation  of  (15).  The  divergence  occurs  in  a  region  that  violates 
this  assumption.  To  remove  the  divergence,  one  should  allow  for  the  possible  mutual 
approach  of  particles  1  and  2  in  (11). 

In  the  classical  gas  of  electrons  considered,  one  never  encounters  the  problem  of 
two  particles  close  together  because  of  the  mutual  repulsions .  This  fact  should  be  borne 
out  in  the  solution  of  the  set  of  equations  (11).  If  discrete  positive  charges  had  been 
included,  real  problems  might  have  occurred  because  of  the  attractive  potential.  The 
occurrence  of  bound  states  and  the  effects  of  very  fast  electrons  would  require  analysis 
that  would  go  far  beyond  the  techniques  employed  here .  In  this  investigation  we  shall 
study  the  effect  of  close  collisions  on  the  two-body  distribution  function  and  shall  thereby 
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remove  divergence  in  the  classical  problem  of  Coulomb  repulsions. 

In  the  past  this  divergence  has  been  handled  by  cutting  off  the  integration  in  (17)  at 
2 

kT/e  ,  the  value  corresponding  to  the  distance  of  average  closest  approach.  Since  the 
dependence  on  the  cutoff  is  logarithmic,  the  final  results  are  not  expected  to  be  much 
in  error  from  a  numerical  point  of  view.  However,  it  is  of  interest  to  see  what  happens 
to  the  distribution  function  for  close  collisions  and  to  see  how  accurate  the  method  of 
cutting  off  may  be . 

In  order  to  proceed  let  us  return  to  the  arguments  that  led  to  the  assignments  of 

orders  of  magnitude  to  the  various  terms  in  (11).  We  found  that  the  requirement  that  the 

pair  interaction  term  be  small  implied  that  the  s  particles  are  mutually  separated  by 

distances  b  »  Let  us  now  imagine  a  set  of  concentric  spheres,  of  possible 

separation  of  twogarticles,  with  radii  e  /kT,  V  e  \jy/kT,  and  The  radii  are  in  the 

constant  ratio  VkTX.j^/e^  =  AMiirn^.  For  a  wide  range  of  plasmas  the  quantity  4TrnX.^ 

is  much  larger  than  one.  For  example,  if  kT  =  100  ev  and  n  =  10^^  cm”^,  then 
3  5 

4irn\jj  ~  10  .  These  spheres  are  quite  distinct,  and  they  are  useful  in  visualizing  the 
process  of  interaction. 

Let  us  label  the  spheres  I,  II,  and  III  in  order  of  increasing  radii.  The  volume  inside 
sphere  I  may  be  considered  forbidden  to  the  particles  because  e  /kT  is  the  distance  of 
closest  approach  for  the  particle  of  average  energy  kT.  Particles  are  allowed  between 
spheres  I  and  II;  but  in  this  region  the  potential  energy  is  larger  than  the  average  kinetic 
energy,  and  thus  the  pair  interaction  term  may  not  be  considered  small.  The  solution 
may  be  considered  correct  from  sphere  II  outward,  and  an  evaluation  of  (17)  shows  that 
the  result  is  exponentially  cut  off  outside  in  agreement  with  the  Debye  theory.  In 
order  to  correctly  handle  the  integral  occurring  in  the  equation  for  Fj,  we  must  consider 
the  possibility  of  particles  occurring  between  spheres  I  and  II. 

Consider  the  number  n(e^X.p,/kT)^^^,  the  probable  number  of  particles  inside 

3  3  -1/2 

sphere  II.  In  terms  of  4Trn\j^  this  number  is  (nkj^)  '  ,  which  is  small.  We  may  argue 
from  a  strictly  probabilistic  point  of  view  that  it  would  be  correspondingly  even  more 
unlikely  that  more  particles  should  be  inside  this  sphere .  We  are  led  to  the  concept  of 
a  "close  collision,"  one  in  which  two  particles  experience  a  short  time  interaction  within 
sphere  II.  In  line  with  the  foregoing  argument,  we  shall  assume  that  the  close  collisions 
are  binary  and  shall  ignore  the  possibility  of  three  particles  occurring  within  this  short 
range . 

The  analogy  with  the  Boltzmann  gas  should  be  mentioned.  For  the  Boltzmann  gas 
we  consider  free  particles  undergoing  binary  collisions.  For  the  Coulomb  case  we  real¬ 
ize  that  the  particles  interact  over  a  long  range,  but  we  use  the  fact  that  the  strong  inter¬ 
actions  occur  only  in  binary  types  of  events.  In  both  cases  those  collisions  that  cause  a 
large  change  in  momentum  are  assumed  to  be  binary. 

We  shall  carry  out  the  solution  to  (11)  under  the  assumption  that  the  interaction  poten¬ 
tial  of  one  pair  of  particles,  1  and  2,  is  not  necessarily  small,  while  all  other  pairs  are 
assumed  to  be  small.  The  Hamiltonian  to  zero  order  for  s  particles,  including  1  and  2, 
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will  be  H°  +  4/,  ,•  Sets  of  particles  not  including  1  and  2  will  be  assumed  to  be  outside 
a  range  corresponding  to  a  close  collision;  their  zero-order  Hamiltonian  will  be  H”, 
and  their  solutions  will  be  assumed  to  be  (14)-{17).  Whereas  earlier  treatments  have 
assumed  a  gas  of  electrons  experiencing  entirely  grazing  collisions,  we  assume  that 
one  pair  really  collides  with  no  limitations .  This  assumption  is  not  as  restricted  as  it 
sounds,  since  we  are  in  reality  saying  that  there  are  many  mutually  separated  pairs  in 
close  collision. 

This  entire  procedure  will  not  get  rid  of  all  divergences,  since  in  the  equation  for 
the  two-body  distribution  function  we  encounter  terms  of  the  form 


I 


)]  dx 


3" 


(19) 


By  the  above-given  procedure  we  correctly  allow  for  approach  of  1  and 

3  —  the  combination  would  entail  a  three-body  collision.  In  section  3.2  it  will  be  seen 
that  this  divergence  can  be  circumvented  in  a  plausible  way.  However,  it  is  reasonable 
to  expect  that,  if  we  kept  this  integral  and  merely  cut  it  off,  the  dependence  on  the  cut¬ 
off  in  the  final  equation  for  would  be  much  weaker  than  the  logarithmic  dependence 
found  in  (17)  because  we  have  carried  the  problem  to  one  more  step  of  accuracy. 

Section  2. 1  will  be  devoted  to  a  discussion  of  the  methods,  operators,  and  notation 
to  be  used  in  solving  these  equations .  Section  2 . 2  will  include,  as  an  example,  some 
discussion  of  the  equations  for  the  large-separation  solution.  In  section  2.3  we  shall 
modify  the  ideas  of  section  2. 2  so  that  they  will  be  applicable  to  the  present  problem. 
Section  111  will  be  a  discussion  of  the  actual  problem,  and  Section  IV  will  be  a  discussion 

of  the  results.  Some  of  the  material  contained  in  section  2. 1  is  taken  from  a  book  by 
1 2 

Bogoliubov  but  is  included  here  since  it  is  not  usually  used . 
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II.  THE  METHOD  OF  APPROACH 


2 . 1  THE  ADIABATIC  HYPOTHESIS  AND  ITS  IMPLICATIONS 


In  order  to  study  the  methods  to  be  used  in  the  solution,  introduce  the  Hamiltonian 
for  the  kinetic  energies  and  mutual  interactions  of  s  particles. 

s 

i<j 

For  simplicity,  we  assume  that  there  is  no  external  field.  The  generalization  is  con¬ 
ceptual  immediately,  since  these  are  one -particle  processes.  However,  the  resulting 
particle  trajectories  are  very  difficult  to  solve.  One  would  not  expect  that  the  evolu¬ 
tions  of  correlations  are  much  affected  by  the  presence  of  weak  external  fields;  thus 
this  model  is  useful  intuitively  for  the  more  difficult  case.  The  strong  field,  in  which 
the  external  field  exerts  more  force  than  most  of  the  interactions,  would  be  complicated. 
But  for  this  simple  case  even  the  most  simple  equation,  the  Vlasov  equation  (4),  is  not 
understood,  since  its  nonlinear  character  becomes  important. 

For  the  whole  system  the  Liouville  equation, 

8D„ 

<20) 


has  as  a  formal  solution 


^N^*l’  ' '  ”  *N*^^  ^  t  ^N^^l*  " "  ’*N*^^*  (21) 

N 

Here,  the  operator  S_^  operates  on  the  particle  coordinates  x^  . . .  and  projects  them 
backward  in  time  t  seconds  on  the  basis  of  the  paths  given  by  their  Hamiltonian;  that 
is,  D^  flows  like  an  ideal  fluid  in  phase  space. 

Note  the  identity  for  an  arbitrary  +5 


'  •  '*N’  ^  [**N'  ^-t***^*!'  • '  '^N'  ®-t  h.  •  •  •  *N' 

N 

This  follows  from  the  definition  of  and  the  fact  that 


at  -t’ 


(22) 


(23) 


Another  property  that  we  shall  need  in  this  investigation  is  the  solution  of  equations 
of  the  form 


^(Xj^,  . .  ,x^,t)  —  [H^;^(X|,  ...  x^,t)J  +  f(x^,  . . .  x^#t)  ■ 


(24) 


=,N 


Let  #(Xj,  . .  'Xj^lt)  e  S_^X(x^,  . . .  Xj^;t);  then,  using  (23)  and  (24),  we  obtain 
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0t  ^  *  ’  '  ’  '  ^2 '  '  *  * 

N 

Multiplying  by  S^  and  integrating,  we  obtain 
X(x^,  . .  .Xj^;t)  -  x(Xj,  . .  .Xjj;0)  =  y 


Xj^;t)  dr. 


(25) 


In  (25)  we  use  the  relation 


eN„N  „N  _ ^ 

^t  ^-t  "  ®t-t  t  -  t„  -  0 
o  o  o 


1, 


:N 


and  consider  it  the  inversion  property  of  the  operators.  Finally,  multiplying  by  S  ^ 
and  resubstituting  we  have 


<j)(Xj^,  .  .  .  x^,t)  -  S_^<t>(Xp  .  .  .  Xj^,0)  +  dr  S_^^__^jf(Xj ,  . .  .  Xj^;t)  . 


(26) 


In  (26)  the  first  term  on  the  right-hand  side  is  the  contribution  of  flow  in  phase  space 
and  the  second  is  the  effect  of  the  source. 

Now  let  us  examine  the  exact  equation  for  F^, 


aFj(x^) 

8t 


=  [H°;F^(x^)]  +  n  ^  [4(j2;F2(Xj,X2)]  dx^. 


(27) 


This  is  obtained  by  integrating  the  Liouville  equation  and  is  equivalent  to  it.  To  use  F2, 
we  must  solve  the  equation  for  F^,  which  involves  knowing  F^,  and  so  forth.  The  advan¬ 
tage  of  this  chain  over  the  original  Liouville  equation  is  that,  if  we  can  break  the  chain 
in  a  physically  sensible  way,  we  can  obtain  a  closed  set  of  equations  and  know  the  pre¬ 
cise  approximations  made  in  departing  from  the  full  Liouville  equation.  The  problem 
is  to  perform  the  break  in  a  manner  that  will  balance  the  physical  and  computational 
reasonability . 

Now  look  at  (27)  in  the  light  of  the  fact  that,  as  mentioned  in  section  1 . 1,  the  var¬ 
ious  forms  of  kinetic  equations  for  F^  have  one  thing  in  common  —  they  can  generally 
be  written 


»F  (X  ) 

“  Aj(Xj,F^).  (28) 


Here,  Aj(x^;F^)  is  a  functional  of  F^ .  The  importance  of  this  is  that  the  entire  time 
dependence  of  the  right-hand  side  of  (28)  lies  inside  F^  and  depends  only  upon  the  current 
value  of  F^ .  This  fact  is  implied  by  having  a  kinetic  equation  for  F^:  that  its  present 
value  is  sufficient  to  predict  its  future. 

If  (27)  is  to  be  of  the  form  of  (28),  then  F^  must  be  such  that 

F2(Xj,X2;t)  =  F2(x^,X2;F^).  (29) 


If  (29)  is  to  hold,  it  must  be  true  for  all  s  >2  that 
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Fg(Xi .  . .  Xg;t)  =  Fg(Xj...Xg;F^). 


(30) 


1 2 

Bogoliubov  has  shown  that  for  a  Boltzmann  gas  this  is  a  very  good  assumption  in 

that  any  initial  F  that  violates  (30)  will  relax  to  the  form  (30)  in  a  collision  time  that  is 
s 

very  short  compared  with  the  characteristic  time  of  change  of  F^ .  Thus  it  is  safe  to 

assume  that  F  is  of  the  form  (30). 
s 

This  whole  argument  breaks  down  for  a  plasma,  particularly  one  that  is  not  spacially 
uniform.  In  this  case  the  collision  time  is  of  the  order  Vt>/v  =‘\^/4irne  =  which 

J-/ 

is  the  characteristic  time  of  change  of  F^  for  a  nonuniform  plasma.  However,  to  get 
a  kinetic  equation  of  the  form  (28),  one  is  forced  to  take  (30)  as  an  assumption  and  to 
look  for  those  solutions  satisfying  this  form  which  will  provide  the  most  general  equa¬ 
tion  of  the  form  (28). 

We  begin  our  investigation  with  a  generalization  of  the  foregoing  equations .  We  shall 
look  for  a  pair  of  equations  which  is  made  up  of  (27)  and  an  equation  of  the  form 

8F,(x,.  X,) 

0t  “  ^2^*1' *2’^1’ ^2^’ 

where  the  time  dependence  of  resides  within  a  functional  dependence  on  F^  and  F^. 
Since  we  know  from  (11)  that 


H°+4-i2;F2 


n  ^  2'*’*J^23'^3^  ^3 


(32) 


we  may  say  in  analogy  with  (29)  that 

and,  therefore,  that 

Fg(Xi,  ...Xg;t)  =  Fg(Xj...Xg;Fj,F2). 


(33) 


(34) 


We  shall  refer  to  these  functional-dependence  assumptions  as  the  "adiabatic  hypoth¬ 
eses"  in  the  following  sense.  For  example,  in  (34),  we  assume  that  F  for  s  >  2  relaxes 

s 

very  rapidly  to  a  form  depending  only  on  the  instantaneous  values  of  Fj  and  F^  •  This 
assumption  is  analogous  to  the  adiabatic  approximation  to  the  time-dependent  perturba¬ 
tion  theory  in  quantum  mechanics . 

This  procedure  serves  the  following  purpose .  We  assume  that  a  kinetic  equation  for 
F2  exists,  (31).  Equation  31  implies  certain  limitations  upon  the  time  dependence  of 
Fj,  (33).  We  shall  find  the  solution  for  F^  which  satisfies  these  limitations  and  use 
this  function  in  the  integral  of  ( 32)  to  obtain  a  general  kinetic  equation  for  F^  •  We  then 
have  a  method  of  investigating  directly  the  possibilities  and  limitations  of  making  another 
restriction  on  F2,  that  is,  that  it  is  of  the  form  F2(Xj,  X2;Fj).  We  shall  make  this 
restriction  and  obtain  a  solution  for  F2  in  this  form .  Since  the  additional  assumption 
that  was  made  to  obtain  (17)  from  (15)  is  (29),  this  procedure  will  yield  a  result  that  is 
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directly  equivalent  to  (17)  except  that  the  divergence  will  no  longer  exist. 

A  look  at  the  equations  to  be  solved  will  show  the  origin  of  the  adiabatic  hypothesis 
and  the  ensuing  statements .  If  the  interaction  ^  is  of  order  g,  a  condition  that  we 
shall  call  Case  I,  we  can  write 


“  S 

-I 

s 

I 

^  “^is+l'^s+l 

_i=l 

(35) 


When  particles  1  and  2  are  allowed  to  collide,  small  and  must  be  included  in 

the  zero-order  Hamiltonian.  In  this  form  we  have  Case  II, 


[«s++12-''s 

+  g 

i 

i<3 

Z  '^is+r^ 

i=l 

where  the  ^  denotes  the  omission  of  the  pair  {l2}. 

We  shall  see  that  if  is  expanded  in  powers  of  g. 


s+1 


dx 


s+1' 


(36) 


(37) 


and  a  perturbation  expansion  is  carried  out,  then  the  zero  solutions  are 
CASE  I: 


F°=rT  F  (x.) 
®  i=l  ^  ^ 


(38a) 


and 


CASE  II: 

7.0 


F“  =  F^Uj.Xj)  n  F,(*,). 


(38b) 


In  the  perturbation  expansions  of  (35)  and  (36),  f“  will  serve  as  the  source  term  for 
1  “ 

Fg  and  we  obtain 

CASE  I: 

4 


ST 


s 

/> 

s 

+ 

2j  *^18+ 1 '^8+1 

i=l 

dx 


s+1 


and 


CASE  II: 
8F^ 


r 

s 

+ 

y  ^^••;F° 

Lj  ^13  S 

Z/  '^is+l'^s+l 

i=l 

dx 


s+1 ' 


(39a) 


(39b) 


With  the  adiabatic  hypothesis  applied  to  (39)  we  can  study  the  time  dependence  of  F  while 
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holding  F°  constant  in  the  source.  If,  upon  solution,  is  found  to  vary  as  rapidly  as 

o  ^  ^ 

F°,  then  the  hypothesis  is  invalid  and  another  method  is  required, 
s 

The  basis  for  the  belief  that  for  a  uniform  plasma  the  hypothesis  is  valid  follows 
from  the  equation  for  F^  in  which  we  expand  and  use  (38) 

8F,  (x.)  r  1  o 

8t  “  [^l'^l^*l^J  ^  "  J  ^^2 

+  g  n  [4-i2'F2<'‘i**2)  dX2  +  .  .  .  .  (40) 

The  first  integral  is  the  effect  of  a  potential  U  whose  source  is  F^, 

V^U  =  -4ime^  J  Fj(x)  dp.  (41) 

For  the  neutral  uniform  plasma,  U  and  the  Poisson  bracket  of  Fj  with  H°  are  zero.  The 

rate  of  change  of  F,  is  of  order  g;  thus  in  equations  like  (39)  we  can  hold  F°  constant 
1  s 

to  the  same  order  in  g. 

In  this  report  we  assume  for  Case  1  and  Case  II  that  the  F  depends  functionally  upon 

o  ^ 

the  functions  appearing  in  F  .  In  section  2 . 2  we  work  out  Case  I  by  assuming  a  depend- 

s 

ence  on  F^ .  This  solution  will  serve  as  an  introduction  to  the  methods  because  this 
problem  is  closely  related  to  the  more  general  one.  Section  2.3  and  Section  III  will  be 
concerned  with  Case  II. 

2.2  THE  LARGE -SEPARATION  SOLUTION 

In  this  section  we  shall  discuss  the  solution  to  those  equations  in  which  all  two- 
particle  interactions  are  assumed  to  be  small,  which  were  referred  to  as  Case  I  in  sec¬ 
tion  2.1.  We  take  (30)  to  be  our  form  of  the  adiabatic  hypothesis. 

To  carry  out  the  solution,  we  must  be  able  to  handle  terms  of  the  form 

it  ^s^*r *2‘  ■  ■ 

for  which  we  know  that  the  derivative  will  operate  only  on  the  F^,  since  that  is  the  only 
place  in  which  t  appears .  For  s  >  2  expand  F^  ip  a  power  series  in  g,  where  g  is  the 
expansion  parameter  described  in  section  1.2,  to  obtain  (37).  In  (30)  the  expansion  of 
F^  will  have  the  effect 

8t  ^  ^  ['^12'^2^*1’ *2'^P'*’^  ^2^*1' *2*^2^^'  '  ‘  *^*2' 

and  we  shall  consider  this  an  expansion  of  OF^/Ot  in  powers  of  g.  For  convenience  of 
notation  we  write  (42) 

- JpL  =  A°(Xj;Fj)  +  gAj(Xj;Fj)+  ...  (43) 
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I 


and  define  the  Aj  by  comparison  with  (43): 

aJ{Xj;Fj)  =  jHj;Fj(x^)  +  n  y  ‘^2' 

a}(Xj;Fj)  =  n  ^  [‘I'l 2*^2 <*!' *2*^1)]  ‘^2* 


(44) 


Let  X(x,.  ...  X  ;F,)  be  any  expression  whose  time  dependence  is  completely  defined 
1  SI 

by  Fj .  For  infinitesimal  variations  in  F^, 

6X(Xj.X2.  . .  Xg'Fj)  =  <|>(Xj,  . . .  Xg;Fj,6Fj), 

where  ^  is  some  new  functional  that  will  be  linear  in  6Fj .  From  this  relationship  and 
(44),  we  have 

"at  x(*i’  •  •  •  *s'^l^  ~  ^ 

=  +(xj....Xg;Fj.A°)  +  g^(xj....Xg;Fj.Aj)  +  . . . 


fx  X  F 


or 


8X  _ 

~st- 


DqX  +  gD^X  + 


(45) 


Here,  D  denotes  an  operator  that  differentiates  with  respect  to  t  (by  operating  on  F.) 

8F  r 

and  then  replaces— with  A^ . 

We  wish  to  apply  these  definitions  to  the  equation  for  F^,  Case  I, 


+  g 

cn  , 

”s 

^'^is+l'^s+l 

Li<j  J 

i=l 

(46) 


In  (46)  expand  SF^/dt  as  in  (45),  and  expand  F^  as  in  (37)  to  obtain 


(D„*gDjtg2...)(F»tgF‘tg2p2,..,)  ,  H°*g^  * ‘''s  ^  • 


KJ 


"I 


^  '^'is+l'^s+1  ■'■^^s+1  ^s+1  •  • 


i=l 


dx 


s+1 ' 


Using  the  first  two  powers  of  g,  we  obtain  the  following  two  equations: 
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(47) 


-I 

■  s 

X  '*'is+l'^s+l 

.i=l 

dx 


s+1 


and 


'  s 

■  s 

D  F^  +  D,F°  =  [h^iF^ 
OS  Is  |_  s  s 

+ 

^  ^is+l'^s+1 

.i<j 

i=l 

dx 


s+1  ■ 


(48) 


Since  in  the  limit  of  infinite  separation  of  all  s  particles  the  motions  of  the  particles 
must  become  uncorrelated,  we  impose  the  following  boundary  conditions: 


— ^  TT  Fi(x.) 

all  -►  00  i=  1 


(49) 


»  0 


i  >  1. 


all  q.  -  q^  -  00 


Here,  goes  to  zero  because  F°  is  assumed  to  have  all  of  the  boundary  contribution. 

By  direct  substitution  and  the  use  of  the  definition  of  the  D  operators  in  terms  of 

r  ^ 

Aj,  we  see  that  the  solution  of  (47)  which  satisfies  (38)  is 


F»(Xi,...x^:Fp.J7  Fi(Xj). 

In  solving  (48)  we  use  the  fact  that,  by  definition  of  D^, 
g  s 

=  °i  n  w  =  1  n 

i=l  J". 


(50) 


s 

“"S  n  Fj(xj)y[*j^^,;F*(x, 

1=1 


^'*S+1'^1^ 


dx 


s+1 ' 


Substitution  of  (51)  in  (48)  yields 


s 

s 

+ 

i  *ijK 

^  “^is+l'^s+l 

i=l 

dx 


s+1 


s 

‘=‘ 


dx 


s+r 


(51) 


(52) 


We  see  that  by  straightforward  substitution  and  application  of  (44)  and  (50),  the  solution 
of  (52)  which  is  consistent  with  the  boundary  condition  (49)  is 
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K— 1 


The  final  equation  for  F2(x^.X2;Fj)  is 


"I  h 


dx. 


3''''^23'^2^*l’ *2'^l^  Fi(x3) 

For  uniform  plasmas,  this  equation  has  been  solved  by  Lenard^*^  and  Dupree. The 
latter's  solution  is  carried  out  in  the  Appendix  by  using  methods  to  be  introduced  in 
Section  III.  This  solution  results  in  the  expression  (17),  which  was  solved  under  the 
adiabatic  hypothesis.  Thus  in  line  with  (18),  the  solution  to  (53)  will  be  denoted 

4<*1'*2'^i'- 

2.3  COLLISION  PROBLEM  ASSUMPTIONS 


In  this  section  we  discuss  the  equations  referred  to  as  Case  II,  (36),  in  which  4'j2 

is  not  assumed  to  be  of  order  g  and  the  expansion  parameter  has  been  inserted  before 

those  terms  taken  to  be  small.  (See  section  1.2.) 

By  selecting  particles  1  and  2  as  the  particular  particles  whose  close  approach  will 

be  allowed,  we  temporarily  destroy  the  interchangeability  of  particles  in  F  .  Thus  when 

s 

we  make  an  adiabatic  hypothesis  for  F  (Eq.  34)  and  take  8F  /8t,  we  must  interpret  the 

s  s 

effect  of  8Fj/8t  and  8F2/8t  differently,  according  to  their  arguments.  Since  in  (36)  we 

assumed  that  all  are  small  for  i  and  j  which  are  not  equal  to  1  and  2,  we  expect 

the  correlations  between  these  pairs  to  be  the  same  as  those  studied  in  section  2.2 

because  the  correlation  of  a  pair  of  particles  is  generated  by  the  past  history  of  their 

mutual  force.  Therefore,  for  these  i  and  j,  we  take  8F,  (x.)/8t  and  8F„(x.,  x.)/8t  to 

be  given  by  the  results  of  section  2.2. 

N  -  s 

In  (10)  we  obtained  the  factor  — ^—before  the  integral  by  summing  over  identical 
particles.  If  in  this  new  interpretation  we  are  careful  to  sum  over  only  the  identical 
ones,  we  will  have  as  an  equation  for 


s+1 


8t  ^  5=  |Hj;Fi(Xi)  +  y  j'[+.g^j;F2(x.,Xg^j)]  dXg 


However,  the  last  two  integrals  have  small  effect  for  two  reasons.  First,  in  the  limit 
of  large  N,  their  contribution  will  become  negligible .  Second,  we  know  that  outside  a 
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range  corresponding  to  a  close  collision  all  F2  are  the  same,  so  that  to  this  extent  they 
can  be  included  in  the  first  integral.  We  have  said  the  or  must  be  small,  just 
as  is  small;  thus  the  contribution  must  be  essentially  the  same.  A  similar  argu¬ 

ment  holds  for  the  evolution  of  F2(Xj,X2). 

This  distinction  may  be  summarized  in  the  following  way.  When  we  write  an  equa¬ 
tion  for  Fg,  the  Hamiltonian  will  contain  ijjj2  io  zero  order  in  g  only  if  the  set  s  con¬ 
tains  both  1  and  2.  We  assume  that,  if  this  zero-order  Hamiltonian  does  not  contain 

the  evolution  of  F  is  the  same  as  that  discussed  in  section  2.2.  This  assumption 
is  connected  with  the  binary  collision  assumption.  We  do  not  complicate  the  picture  by 
including  close  collisions  of  mutually  separated  pairs  and  we  explicitly  exclude  3 -body 
collisions . 

With  this  in  mind  we  are  able  to  interpret  differentiations  of  the  form-^F  (x,,  .  , .  x  ; 

0t  s  1  s 

F^,  F2).  When  the  derivative  operates  on  Fj(x^)  or  F2(Xj^,  x^),  we  apply  defined  in 
section  2.2.  When  the  derivative  operates  on  F2(Xj,X2)  we  must  define  new  A2.  By 
comparing  (31)  with  the  equation  for  F2(Xj,X2), 

®2''‘'*'l2'^2^*l'*2^]  ”  y  ['J'i3++23'^3**1'*2’*3'^1' ^2^^  *^3' 

expanding  in  (54) 

F3  =  F°+  gF^  +  g^F^  +  ... 


2'  1'  2 
at 


and  expanding  in  (31) 

A2<*i«X2'^1'^2)  =  A2(Xj,X2;Fj.  F2)  +  g aJ(Xj. X2;Fj.  F2)  +  g^  . 
we  obtain 

*2'^r  ^2^  =  |H2++j2iF2(Xj,  X2) 

+  n  j*  [^j3+^;23;F°(x3.X2,X3;Fj,  F2)]  dx3 
^2^*1' *i'^r  ^2^  =  n  J  jq/j3-H|<23;F3(Xj,X2,X3;Fj,  F2)j  dXj 


(55) 
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III.  COLLISION  PROBLEM  SOLUTION 


In  section  3.1  we  shall  use  the  results  of  section  2.3  in  an  equation  for 
F^.  In  section  3.2  we  shall  solve  this  equation  and  use  it  in  a  solution  of  the 
equation  for  F2 . 

3.1  THE  EQUATION  FOR  F^ 


In  order  to  derive  an  equation  for  F2(x^,  x^,  x^iF^,  F^),  we  start  with  Eq.  36  in  which 
we  wish  to  obtain  an  expansion  for  F^  which  is  analogous  to  (37).  In  (36)  we  no  longer 
require  that  i)'j2  order  g,  and  we  assume  that  its  time  dependence  occurs  through 

a  functional  dependence  on  Fj  and  F2.  Expand  8/dt  as  in  (45),  F^  as  in  (37),  and 
then 


(F^g  F‘t. . .)  .  F‘t. . . 

+  g 

~S  1 

r 


*^3 


l 

i=l 


'^iS+l'^S+1 


S  +  1 


dx 


s+1  ■ 


Equating  the  first  two  powers  of  g,  we  obtain 


-5 

y  4-  ;f° 

/j  ^is+r  s+1 

i=l 

dx 


s+1 


and 


■  s 

*  s 

K+'I'i2-f‘_ 

+ 

Z  +irF“ 

*^is+l'^s+l 

_i=l 

dx 


S+1  ■ 


(56) 


(57) 


Since  in  the  limit  of  infinite  separation  of  all  s  particles  except  1  and  2  the  motions  must 
become  uncorrelated,  we  impose  the  following  boundary  conditions 


F-(Xi,...x^;Fi.F2) 


Fg(x^,  .  .  F2) 


- ^ ^  F,(x,,x,)TT  F,(x,.) 

allqj^i-oo  ^  ‘  i=3  ‘  " 


except  12 


(58) 


i  >  1. 


all  q^  -  q^  -  00 
except  12 


The  boundary  condition  for  F2(Xj,X2)  will  not  be  introduced  here  because  it  is  not  needed, 
since  F°  has  F  (x, ,  x„)  in  it.  The  boundary  conditions  are  used  only  to  show  the 

S  6  i  b  ^ 

asymptotic  form  of  F^. 

The  solution  to  (56)  which  matches  these  boundary  conditions  is 
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Fj)  .  F2(*i.*2)  n  Pi(Xi). 


Here,  we  have  used  the  facts  that 


and 


*2^  =  A2(Xj.X2;Fj,  F2) 


DoFi(Xi)  =  Ai(x.;Fi). 


To  solve  (57),  we  must  know  the  function  D,F  .  Using  the  section  2.2  prescription 

1  s 

for  Dj,  we  have 

s 

n  D,F2(x,,x,)  t  F2(x,.X2)  1  ft  Fj(*.)  DjFjlXj) 


.  1=3 


"U  ^  [‘*'ls+l''’‘*'2s+l'^3^*l’*2'^s+l’^l' ^2!] 

s 

+  nF2(x^,X2)  ^  n  y  ['<-is+l4<V*s+l'^l) 


1=3 
i=3  •*  r 


s+1 

s+r 

•  th 


(59) 


The  implications  of  the  discussion  in  section  2. 3  are  used.  Since  the  i  particle  is 

assumed  not  to  interact  closely  with  the  other  particles,  the  function  (5^) 

is  the  large -separation  function  of  section  2.3. 

Using  (59)  in  (57)  we  arrive  at  the  following  expression  for  D  F^: 

o  s 


*  s 
v~ 

Z  • 

..Xg;Fj,F2) 

I  ‘»'is+i'^Li<*r--Vi'^i'^2> 

-i=l 

~  "  i^3  ['*'ls+l''''*'2s+l'^3^*l'*2’*s+l'^l' ^2^ 


dx 


s+1 


S  g 

-nF2(Xj,X2)  ^  n  Fj(x.)  j[g;.g^j;4(x., 

'ii 


^s+P 


*^s+r 


By  direct  substitution  the  solution  to  (60)  is 
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*2)  I  ft  Fj(x^)  4(Xj.  X.) 


a  s 

^  I  JI  F^(Xi.X2.x.;Fj,F2). 


i=3 


s 

A 

#i 


,1  ... 


Here,  F^  is  the  solution  to 


^0^3^*!'  *2'  *3'^r  ^2^  [^3‘'''*'l2’^3^^1'  *2’  ^2^_ 

+  n  j*[4<i4+4'24+'<'34'^l<*4)Fi<*l'*2'^3'^l'^2)]  *^4 
■'■  ^  J['*'34'^l  (x2)F3(Xj,X2,x^;Fj,  F2)j  dx^ 

"  I  ['i'l4+'*'24'^2<'"l'*2)4<*3''‘4'^l) 


dx . . 
4 


(61) 


The  first  integral  is  zero  for  a  uniform  plasma,  since  it  is  the  effect  of  the  potential  U 
in  (41). 

For  convenience  we  put  (6I)  into  another  form  and  thus  give  an  analytical  meaning 
to  D^. 

o  00 

Since  for  the  uniform  plasma  the  integral  terms  in  A2(Xj,X2)  and  Aj(Xj)  are  zero, 

we  have 


^o^2^*r*2^  "  [^2^'^12'^2^*1' *2^j 


(62) 


Now  introduce  the  operator  S_^  that  projects  particle  3  backward  t  seconds  along 
a  path  given  by  the  free-particle  Hamiltonian.  For  any  function  ^(x^). 


-|;Si,+(X3)  =[h°;s1^+(X3|  . 


(63) 


Likewise  introduce  the  operator  S_^  that  projects  particles  1  and  2  backward  along  the 
paths  given  by  their  Hamiltonian  including  interactions .  In  operating  on  the  coordinates 
of  a  colliding  pair  at  x^  and  x^,  this  operator  will  produce  the  coordinates  Xj(-t)  and 
X2(-t).  In  analogy  with  (63),  we  have 


S_^+(x^,  x^)  =  H2;S_^^(Xj,  X2)  . 


(64) 
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Using  FjUj)  and  F2{x^,X2)  in  (63)  and  (64),  we  have 

f  si^Fi(x3)  =  [Hj;si^Fi(x3) 

87  =  |H2+4'j2iS_^F2(Xj,X2)  . 


(65) 


Comparing  (62)  and  (65),  we  have  the  identity 


(66) 


1  2 

that  is,  if  in  the  functional  dependence  we  use  S_^Fj  and  S_^F2  instead  of  Fj  and  F2, 
the  operator  can  be  replaced  with  8/8t  because  the  effect  is  the  same.  By  using  this 
relationship,  (6I)  becomes 

87  ^3^*1'*2'*3'®-t^1' “  [*^3‘‘’'*’l2'^3  ^*1’ ^^2' ^3'®-t^1' ®-t^2)] 


+  n  y  [+34<F^(x3)F^  t^,X2,x^<F^,  s2^F2)]  dx, 

+  |»|<j3+»l'23iS_^F2(Xj,X2)S_^Fj(x3)j 


+  n  j'[4'i4+4<24;sf^F2(Xj,X2)£j(x3,X4;si^F^)]  dx^. 

(67) 

1  2 

where  we  have  replaced  Fj^  by  and  F2  by  S_^F2. 

Note  that  in  (67)  t  is  not  the  time  variable  t.  Here,  t  occurs  as  a  parameter  inside 

F,  and  F,,  and  t  is  a  dummy  variable  introduced  to  give  analytical  meaning  to  D  . 
lb  o 

Equation  67  holds  for  any  value  of  t;  we  shall  pick  the  value  of  t  that  is  most  convenient 
to  us .  Because  of  the  adiabatic  hypothesis  and  the  implications  of  it,  t  lost  its  position 
as  a  variable  and  became  a  parameter  in  the  equation  for  F^.  The  t-dependence  has 
become  a  functional  dependence  on  F^  and  F2. 

In  order  to  make  the  functional  substitution  clear,  (67)  is  more  general  than 
is  necessary.  Since  we  have  assiuned  that  F^  has  no  spacial  dependence,  we 
can  use 


S^^F^(?3)  =  F^(|?3)  (68) 

throughout  this  report. 

3.2  SOLUTION  OF  THE  EQUATION  FOR  F^ 

We  shall  now  solve  the  equation  for  F3  for  use  in  the  equation  for  F2  without 
expanding  F,  in  g;  however,  we  must  expand  F-  in  order  to  obtain  its  solution.  Some 
of  these  operator  techniques  were  developed  by  Dupree . 
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To  proceed,  go  back  to  (79)  and  use  (80),  to  obtain 
^F3(xj,X2,X2;Fj,  S_^F2)  =  j^H2+4(j2*F3(xi.  *2' ^3'^!' 

+  n  [^34;Fj(p3)Fj(xj,X2,X4;Fj,sf^F2)]  dx^ 

+  |4'j^3+»l'23iS_^F2^*l'*2^Fi(?3)J 

+  n  j'[+i4+4'24'S!TF2(Xj,X2)4(x3,X4;Fj)]  dx^. 

For  the  moment  abbreviate 

F3^*1'*2'*3'Fi' ^-T^z)  = 

where  the  other  dependences  are  understood.  Using  (22),  we  may  then  write  (69) 
^sJ,F‘(r)  t  Ls2,f‘(t)  = 


(69) 


(70) 


Here,  the  operator  L  acting  on  any  function  +(x3)  is 

^(Xj)  -[h?;+(X3)]  J  4.3^+(x,,  dx^,  (71) 

->  8p3  8q3 

2 

and  +(t)  is  the  source  term  made  up  of  the  last  two  terms  of  (69).  Note  that  and  L 
commute  since  they  operate  on  different  coordinates.  We  call  L  the  Landau  operator 
since 


8f(x3,  t) 


8t 


+  Lf(x3,  t)  =  0 


is  the  equation  involving  L,  which  has  been  discussed  in  detail  by  Landau . 

Writing  out  the  arguments,  we  obtain  the  formal  solution  of  (70): 

FjCi.X2.X3;Fj,sf^F2)  =  e"^""  dr'  e^'"'  X2. X3;Fj.  sf^,F2) 


+  e"^"^  sf^4»(Xj,X2,X3;Fj,  F2). 


(72) 


where  we  use  the  fact  that  L  does  not  depend  on  t  when  F,  has  no  'q-dependence. 
2  ^ 

In  (72)  let  F2  S^^F2,  use  the  inversion  property  of  the  S_^  operator,  and  let 
t'  •  T  -  t',  to  obtain 


F3(*i«*2'*3*Fi'^2)  =  r  e  sl^i*(xyX^,XyF^.S^,F^  dr' 


+  sf^+(Xj,X2.X3;Fj,  F2). 


(73) 
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The  left-hand  side  of  (73)  is  independent  of  t,  and  thus  we  are  free  to  pick  t  arbitrarily. 
We  pick  T  =  00  and  thus  remove  the  initial  condition  term  because  by  (58)  must  go  to 
zero  for  infinite  separation  of  1  and  2. 

We  then  have  a  new  understanding  of  the  adiabatic  hypothesis  and  the  expansion. 
We  shall  see  that  (73)  implies  that,  to  this  order  of  g,  the  correlations  are  calculated 
by  integrating  along  unperturbed  (that  is.  zero-order)  orbits  while  and  F^  are  held 
constant.  If  F^  and  F^  are  known  to  change  in  times  that  are  comparable  to  the  time  of 
build-up  of  correlations,  for  example  in  a  nonuniform  plasma,  then  this  analysis  is 
incorrect.  The  dummy  variable  t  gives  us  a  way  of  studying  this  mathematically. 

We  must  interpret  the  meaning  of  the  operator  e”^^  in  (73).  For  any  ♦(x^),  e~^^ 
satisfies 

e  <|)(X3)  +  L  e"^"^  ^{x^)  h  o.  (74) 

For  the  moment  let 


h(x2,  t)  s  e  <t)(x2) 

for  any  <t>(x2),  and  let  us  represent  by  h  the  Fourier-Laplace  transform 
h(P3,  k,  9)  =  dr  e~^'^  ^(pj.  k,  t) 


(75) 


where 


-  -  C 

Mpj.  k,  t)  =  J  dq^  e  ^  Mp^.q^.r). 

A 

As  shown  by  Landau,  (75)  reduces  (74)  to 


2. 

mu  1 

^  P  ^ 

(<r-ik.v,)  h  + - 5— k  • 


8F  (p  )  p.  ^  ^ 

— 1  \  h  dp'  =  h(p  k.  0) 

8p,  J  3  3 


(76) 


or 


I 


dP3  = 


h(p3,  k,  0) 
or  -  ik  •  V. 


dp. 


1  +  m  ■ 


2. 
u  1 
P 


I 


BP' 


■dp' 


<r  -  ik  . 


(77) 


(We  are  actually  interested  only  in  the  integral  of  h  over  dp3.)  Since  Re  a  >  0  in  the 
definition  of  h  in  (75),  the  integral  in  (77)  is  defined.  To  stay  above  the  point  (-L<r), 
we  take  the  same  integral,  with  the  P3 -integration  deformed,  to  be  the  analytic  continua¬ 
tion  of  the  function  Re  (r  <  0.  All  of  this  computation  is  identical  to  Landau's  work. 
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Finally,  inverting  the  transform,  we  obtain 


1 

(2ir)^  i 


'ooi+p 


pooi+ 

\  d<r  e 

*^-ooi+p 


dk  e  ^ 


h(p,,  k,  0) 

- 2__dp 

0-  -  ik  •  Vj 

1  +  L  +  (o-) 


3 


where  we  define 


(78) 


L^(a)  H  m 


^  8F  (p'  ) 
k - 


(79) 


The  plus  and  minus  signs  indicate  that  the  contour  is  to  pass  above  or  below  the  singU' 
larity . 

Resorting  to  the  definition  of  h,  we  have 


I 


e  <t)(xj)  dffj 


(zttr  i  *^-1 


pioo+p 


do-  e 


•ioo+p 


I  ^ 


--  I 


<j>(k,p  ) 

dp. 


-ik-q^  0-  -  ik  •  v^ 


1  -t  L^(o-) 


This  implies  that 


J  dPj  F3(Xj,X2,X3;Fj,  F2) 

^  1  dr  e  <|)(s_^Xj,  X2.Xj;Fj,  F2) 


(2v)^  i  ^0 


dr  \  dor  e  \  dk  e 
d-goi+P 


I 


-ik-a. 


*  .  o-T  r  ^  0-  -ik  •  V, 


1  +  L^(<r) 


(80) 


where 


4>(s_^Xj,  X2,  k,  PsxFj,  F2)  ^  dq^  e  +(s_^Xj,  X2,  x^.F^,  F2)- 

In  (80)  we  have  used  the  fact  that 

X2,  X3;F^,  S2F2)  =  X2,  X3;Fj,  F2). 


(81) 


2 

This  follows,  since  S_  is  defined  to  operate  inside  the  function  wherever  x,  and  x, 

^  2>  X  b 

occur.  In  (81)  we  interpret  the  S_^  as  operating  on  those  x^  and  X2  lying  outside  the 
function  F2(x^,  X2)  • 
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To  continue  we  must  insert  the  form  of  Remember  that  we  defined 


Using  these  definitions  we  obtain 

x^j  kj  ^2^ 


(82) 


8F, 


'*13^2^*!'  *2^  ^3^3^  ^  ” 


.(Xi^x^)  p  , 


®Pi 


I 


f^(?3.-p4,  k)  dp^  )  e 


ik-q, 


In  (83)  the  notation  (l'^^2)  means  that  the  previous  term  is  repeated  with  1  and  2  inter¬ 
changed. 


^  g_8 _ 8_ 

8p.  8p^ 

and 

^2^3' ^4'^^  =  ^  dq  4(53,  P4,q)- 

If  qj(-T)  is  the  position  of  particle  1  at  -t,  given  that  it  was  at  q ^  at  t  =  0,  with  the  mech¬ 
anism  of  the  motion  governed  by  the  two-body  collision,  and  similarly  for  q2(-T),  then 
(83)  becomes 


dj3F2(Xj.X2)  FjCpj)  +  n— i-j  4(P3.p4,E) dp^j  e  ^  +(hM-2)j., 


(84) 


Let  us  now  discuss  the  interpretation  of  some  of  the  integrals  appearing  in  (80).  Note 
that  the  T-dependence  in  ^  is  in  the  form  of  an  exponential  with  an  imaginary  argument. 
Thus  the  r-integration  can  proceed  only  if  p  is  negative.  Since  the  inversion  in  the 
(T -plane  must  go  to  the  right  of  all  singularities,  p  can  be  negative  only  if  all  of  the  zeros 

of  the  denominator  lie  to  the  left  of  the  imaginary  axis.  For  a  wide  range  of  F,  (p^)  this 

I  I  13  ^ 

is  true  as  long  as  |  k  |  >  0 .  Thus  we  exclude  an  infinitesimal  region  from  the  origin  of 

k.  Under  these  restrictions  we  may  proceed  with  the  x -integration,  followed  by  the  o-, 

then  by  the  k. 

Those  distribution  functions  that  give  zeros  of  1  +  L^((r)  in  the  right -half  plane  pre¬ 
sent  a  difficult  problem.  For  these  situations  the  past  is  not  damped  out,  and  the  cur¬ 
rent  value  of  F3  must  depend  on  the  initial  conditions  —  a  situation  that  is  incompatible 
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with  the  adiabatic  hypothesis. It  would  appear  that  the  hypothesis  is  invalid  for  plas¬ 
mas  containing  these  instabilities.  In  order  to  proceed  we  shall  limit  ourselves  to  the 

1 3 

more  well  behaved  functions  .  As  Baccus  shows,  these  are  all  single-humped 
momentum  distributions. 

A  physical  understanding  of  why  some  distributions  F^(p)  give  rise  to  damped  solu¬ 
tions  and  others  do  not  can  be  seen  from  the  arguments  used  to  explain  Landau  damping 
and  plasma  instabilities.^^  Because  the  shielding  integral  (the  L  operator)  in  (70)  is  the 
same  as  that  considered  by  Landau  in  his  use  of  the  Vlasov  equation,  we  obtained 
Landau's  characteristic  denominator  1  +  L^((r)  and  can  use  his  analysis  of  the  study  of 
waves  to  study  the  growth  of  correlations. 

For  the  partial  wave  in  the  analysis  of  wave  motion  in  a  plasma,  those  particles 
traveling  at  the  phase  velocity  of  the  wave  will  see  a  constant  field.  If 


k 


(p) 

8P  ^ 


v=<r/ik 


is  negative,  more  particles  will  be  speeded  up  than  are  slowed  down,  and  damping  of 
the  wave  will  result.  If  this  derivative  is  positive,  energy  will  be  fed  from  the  particles 
to  the  wave,  and  a  growing  wave  will  result. 

In  the  problem  considered,  the  plasma  is  homogeneous  and  the  k-analysis  refers  to 
the  coordinates  and  ^2  ~  ^3‘  waves  we  had  damped  or  increasing 

energy,  in  the  problem  considered  we  have  decreasing  or  increasing  correlation. 

Using  the  form  of  4*  in  (69)  and  the  meaning  of  the  S_^  operators,  we  have 


♦(s_^Xj,  X2.  Xj,  Ff,  F2)  =  |^S_^(i|<j3+»|'23)'^2^*1' 

+  n  J[s^^(+14++24>'^2<*1'  *4^ 


dx^. 


(85) 


Throughout  this  analysis  we  have  assumed  that  we  never  have  a  three -body  collision, 
that  is,  that  the  interactions  of  all  pairs  except  1  and  2  are  of  order  g.  This  assump¬ 
tion  implies  that  the  distances  Iq^-q^ |,  |qp=qj|.  jq^-q^l,  and  |5p^4l  are  large  in  the 
sense  of  the  g  expansion.  The  S_^  operators  in  (85)  have  the  effect  of 

"^3)-  (86) 


The  change  in  will  become  important  only  when  |qj(-T)-q^  |  becomes  of  the  order 
of  157^3!  ■  assumption  this  change  takes  a  long  time,  a  time  that  is  sufficient  for 
1  and  2  to  be  far  apart.  For  example,  we  assume  that  (57^3!  157^3 1 

At  this  separation  the  potential  2  magnitude 
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which  is  much  smaller  than  the  average  kinetic  energy  kT.  Thus  for  those  values  of  t 
in  (86)  which  have  any  effect,  1  and  2  can  be  considered  to  be  interacting  only  weakly. 

In  (86)  we  will  specifically  use  the  fact  that  as  is  projected  backward  in  time,  the 
difference  |  is  relatively  insensitive  to  the  change  in  q^  while  q^  is  in  the  imme¬ 

diate  vicinity  of  ^2-  A  similar  argument  holds  for  4^2  3*^7  interchanging  1  and  2. 

These  statements  can  be  made  quantitative  by  introducing  in  (86)  the  operator 
that  projects  1  and  2  backward  with  zero  interaction,  that  is,  along  straight  lines  and 
with  constant  velocity. 


s!,+i3(qr^3l  ■  ^4,,  3(57-53)  =  S^,5^+,3(q-v7T=53). 


(87) 


As  we  shall  see  below,  for  large  values  of  t,  the  product  becomes  stationary, 

independent  of  t.  This  fact,  coupled  with  the  preceding  discussion,  will  be  used  to 
approximate  (87). 

To  understand  the  foregoing  assertions,  consider  the  operation  on  the  velocity  of 
1  or  2 


lim 
T  -•  00 


-T 


V.  . 
1 


This  will  approach  a  constant  limit  defined  as 


lim 
T  ••  00 


V. 
-T  1 


=  Vi(Xi,X2) 


By  its  definition  V^(P^=mV^)  is  the  velocity  (momentum)  that  the  i  "  particle  had  in  the 
distant  past  before  undergoing  the  two-body  collision,  given  that  the  colliding  particles 
have  current  coordinates  and  X2*  For  this  reason  we  shall  refer  to  it  as  the  precol¬ 
lision  velocity  (momentum) . 


Equation  88  proceeds  at  such  a  rate  that 


lim 
T  00 


-0. 


This  condition  ensures  the  existence  of  the  limit 


where 

5.(Xi,  X2)  ^  (^i'S-T^i) 

To  find  the  position  Q^(Xj,X2),  start  at  the  point  q^,  project  backward  t  seconds  along 
the  trajectories  of  a  two-body  interaction,  and  project  forward  t  seconds  with  the  con¬ 
stant  precollision  velocity  V. .  In  the  limit  t  00,  this  procedure  will  give  the  position 
that  the  i particle  would  have  had  with  no  collision  interaction.  For  this  reason  we 
call  the  undeflected  position. 
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Now  consider  for  any  <t){Xj.X2) 


'  ♦(si,Pi.  t  (s!,v,-sf^,v,)  dr',  q^  t  (si^T^-si,,,,)  dr^ 


.♦(Pj.Pj.Qj.Qj) 

=  <t>(Xj,X2). 


We  define  the  6 -dimensional  vector  =  {P^,  Q^}. 

Notice  that  the  transformation  from  (x^,  X2)  to  (X^,X2)  is  a  contact  transformation, 
since  it  is  obtained  from  a  product  of  transformations  governed  by  Hamiltonians.  This 
property  will  be  useful  because  of  the  invariance  of  the  Poisson  brackets  with  respect 
to  a  contact  transformation. 

2  *2 

To  approximate  the  product  S_  S  on  any  function  ^(x,,  x^),  we  may  expand  the  result 

T  X  1  fa 

around  t  =  00: 


S^<t.(x, ,  X,)  = 
^  T^'  1'  2* 


•2 

*2^  +  higher-order  terms. 


Application  of  this  to  4'i3(Qj“Vj^T-q^)  yields 

s!,s;+„(q7v^3)  S.|-j3(Q,-V,t^3) 

t  (S^q,-(VV,r))  •  ^  4.i3(«,-V,t^3)  t  . . .  . 


The  coefficient  of  the  second  term  is,  as  a  function  of  r,  the  deflection  of  particle  1 
along  the  path  of  its  collision  with  2.  By  using  the  specific  form  of  t|i.  the  ratio  of  the 
second  to  the  first  term  is  of  order 


|s!,?i-(Qi-Vit||  |q^_ 


(89) 


The  quantity  |Qi“‘13I  is  large  for  the  same  reason  that  jq^-q^  (  is  large,  and  it  is  of 
the  same  order  of  magnitude.  The  ratio  is  then  small  for  all  situations  except  those  in 
which  1  and  2  are  widely  separated  but  have  had  a  collision  in  the  past  (that  is,  (q^-Q^  | 
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is  large  enough  to  be  in  the  range  [Q^-q^l).  For  1  and  2  close  together  and  1  and  2 
experiencing  a  close  collision,  the  ratio  (89)  is  small. 

We  shall  ignore  the  exception  noted  above  and  assume  that  the  ratio  is  small  for  all 
situations  of  interest.  We  shall  see  that  those  situations  for  which  (89)  is  not  small 
have  negligible  contribution  to  the  problems  of  our  interest. 

These  approximations  are  in  the  source  term  of  F^.  If  we  were  to  consider  correc¬ 
tions,  they  would  occur  as  new  (small)  source  terms  and  be  additive. 

We  shall  then  use  the  fact  that  the  approximation 


S^x't'i  3<^3)  =  +1  3<Qi  -^1  -^-53)  (90) 

2 

is  a  valid  outcome  of  the  g  expansion;  and  similarly  for  Since  in  the  integra¬ 
tion  over  dx^  in  the  source  term  ()>  we  assumed  that  jq^-q^  |  and  (  are  large  (that 

is,  we  cut  off  the  integral),  we  can  make  similar  expansions  for  and 

This  step  represents  the  real  departure  from  the  analysis  when  all  particle  inter¬ 
actions  are  assumed  small.  Rather  than  the  approximation  (90),  these  treatments  have 
implicitly  used  the  approximation 


=1 

— T  — T 


or 


=4-i3(qi-ViT^3). 


(91) 


By  the  discussion  above  we  see  that  these  approximations  are  valid  only  when 
weak  at  present  and  in  the  past.  Use  of  (91)  and  the  following  procedure  would  lead 
directly  to  the  large -separation  solution. 

We  may  then  use  these  approximations  in  (84).  Since  we  effectively  make  the  trans¬ 
formation  q^  —  Qj  inside  »|<,  to  retain  the  form  of  the  Poisson  brackets  we  must  trems- 
form  d/dp^  ••  8/8P^.  With  these  changes,  (84)  becomes 


♦^®-T*l'*2'^'P3'^r  ^2)  "  j^2  ^  ■■^^°13^2^*1'*2^  ^1^P3^ 

8F,(x..x,)  f  _.  \  ilc-(Ci-^.T)  1 

+  "-^-  -0  ^2<p3'  P4' 


(92) 


where  D,  ,  s  - ^  and  _  s  - 

8P^  8?j  8P2  8?3 


Using  (92)  in  (80),  we  obtain  another  expression  for  F,, 


1  r_^3_ 

1  +  L  .((t)  O'  -  ilc't 

T 
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^  r/*  r-l  -  \  'I 

^1<P3)  "  ^  -  3  I  ^  +(1*-2)|. 

(93) 

If  p  <  0,  we  may  carry  out  the  t -integration,  which  will  bring  down  (cr-ik-Vj)  ^  and 
((r-ik-^2)  ^  two  terms.  Note  that  the  approximation  (90)  was  made  to  make  the 

T-integration  possible.  For  those  functions  Fj(p2)  for  which  1  +  L^((t')  has  no  poles  in 
the  right-half  plane,  the  <r -integration  can  be  closed  in  the  right-half  plane  and  will  only 
enclose  poles  at  ik-V^  and  ik.V^.  Carrying  out  these  two  steps,  we  obtain 

J  F2(Xj,X2,X3;Fj,  F2)  dp^ 


^  -ik-q. 


+  k-(Vj-V3) 


^13^2^*r*2^  Fi(p3) 


8F  (X  x  )  p  ^  ^ 

+  n  —  \  f2(p,.'P4.  k)  dp  - 2 +  (l-«*2)  L  (94) 

8Pj  J  i  4  Vl+L^(ik-Vj)  J 

The  ^  means  that  the  p^-integration  is  to  stay  above  the  poles,  since  (-io")  was  to  be 
below  the  axis. 

We  wish  to  use  this  result  in  an  integral  of  the  form 

n  y  j+i3(?7^3^‘'''^23^V^3^’*^3^*l'*2'*3'^l'  ^2^  *^3' 

where  we  are  faced  with  the  integration  of  a  coordinate  that  is  assumed  to  be  far  removed 
from  and  Accordingly,  we  may  approximate  under  the  integral 

y  ['*'1 3^^1  ■^3^'^'*'23^‘^2“^3^'^3^^1 '  *2'  *3'^1 '  ^2^]  ‘^3’ 

This  approximation  can  be  analyzed  exactly  as  before,  that  is, 

‘‘'13<^3)  =  '''13<^3>  +  <^1>  '  +  •••  • 

8qi 

The  error  term  is  of  order 
1^1 1 


i'll -^31 


which  has  the  same  characteristics  as  (89)  ■ 

Putting  (94)  into  (95)  yields  the  interaction  integral  in  the  equation  for  F^: 
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dF^ix^.x^) 

8t 


(-!: - ‘  J  ^  (di3I'2I>‘i.>‘2I  flip 

tap^  1  +  L^(ik-V^)  ‘^+  k.(Vj-V3)  \  13  2  1  Z  1 

ep^lXj.x^)  n  ,  ^  r,  ,^\ 

"" - 8?; - j  2<P3'  P4'  ‘^P4/ 

1 _ r  _  ^ 

ik-V,))  k- 


ik.(Q2-Q3) 


8P2  (l+L^Uk  Vj))  ‘^+  k-(V^-V3) 

8F2(Xi.X2)  P  _i  ^  ^  ^  ^  \1  I 

n - - ]  f2(P3.P4.k)  dp^jl  +  (l-*2)j. 


^  -  k  •  6^3F2(Xj.X2)  F^(P3) 


(96) 


This  is  an  equation  for  8F2/8t  to  first  order  in  g.  It  is  complicated  but  will  yield 
some  interesting  insights.  The  integral  terms  can  be  simplified  somewhat  by  examining 
the  equation  for  fi . 

^  11 
Let  us  rewrite  (53)  for  a  uniform  plasma  by  using  F2(x3,  Xj;Fj)  —  f2(x3,  x^)  and  the 

fact  that  f2  depends  only  on  the  difference  of  spacial  coordinates  q  =  q^^=^3 ;  we  obtain 


1  (P3  Pj)  ft  1  ^  ^ 

Do4(P3,Pl.q;Fj)+^j5-i-.^4(P3.Pl,q;Fi) 

8q 


8F,(p,)  g  p  ,  ,  ,  , 

+  n  ^  ^  \  dMr+(g=qS  f2(P3.^=.q  ;f.) 

8p^  8q  J  5  1 

®F  (p  )  8  p  j*l,, _ ..L  ,1,-.  ...  -.L 

-  n - - - \  dp  dq  ij>(q-q  )  fo(P,..  P, ,  q  ) 

8p,  8q  5  251 


=  4: +(5)  •  d3iF^(^3)  Fj(p3).  (97) 

8q 

Again  we  introduce  fi(p,, p, , q;S^  F  )  and  use  the  fact  that  D  -►  8/8t.  We  are  again 
free  to  pick  the  value  of  t,  since  it  is  a  dummy  variable.  However,  since  t  appears  in 
the  unknown  function,  its  value  must  be  selected  carefully.  We  shall  pick  t  -*  00  but  do 
so  in  the  Laplace -transform  space  by  letting  the  transform  variable  o-  —  0.  This  will 
also  set  equal  to  zero,  as  it  must  since  D^F^  =  0.  This  technique  can  be  used 

because  of  the  fact  that  for 
poo 

g(o’)  =  \  g(T)  e  dr,  (98) 

•Iq 
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and 


The  initial  condition  term  will  also  disappear  in  the  limit  o-  —  0. 

Putting  (98)  into  (97)  and  carrying  out  a  Fourier  transform  in  space  as 


we  have 
lim 


4,rne2ii^  (p, )  p 

- dP54(?3-?5'") 


.i  2  . 

4we  1 


—  I'  -  d31^1<P3>  ^1<^)- 
Dividing  by  (€-ik-  (v^-v^))  and  then  integrating  over  p  yields 


.  8Fj(?3) 


J  fl(P3,p,,k)  dp3  =  d?3  J 


aPi  k  .  +  ic  ^  J  f 


Rearranging  terms  and  taking  the  limit  e  -  0,  we  have 
2 


”  j*  ^2fP3'p6-^>  ^5  jdPj 

1  +  L^(ik.Vj) 

We  may  ua.  (100),  together  with  v,  -  V,  aodt^  -  7,,  i„  |,6)  to  obtaio 

VV]  ■ .  JJ5.4  [(.  f  ,;e,.  t„  h  W 


in  (75). 


(99) 


(100) 
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Fj(Pj)  L^(ik-Vj)^ 


1  +  L^(ik-Vj)  J  ^  8F^(P^) 


1  dF^ix^.x^) 

1C 


8P, 


8P. 


F2(x^,X2)  L^(ik-Vj) 
1  +  L^(ik-V^) 

F^(P^)  L^Uk>V^)l 

1  +  L^(ik-Vj)  J  ^ 


+  e 


ik  ■(Q2-Qj) 


[n  ^  f2(P3,  Pj.iO  dp. 


- - =-  fk  •  F  (x  .X2)\ 

8F3(Pj)  8P2  \  8P^  ^  V 


8P, 


L^(ik-Vj)  8F2(Xj.X2) 


l  +  L^(ikV^)  8P2 


+  U« 


2) 


(101) 


This  form  of  the  equation  for  F2  is  still  very  complex;  however,  with  it  we  can 
easily  obtain  results  to  be  compared  with  (17).  Obtaining  these  results  will  entail  making 
an  adiabatic  hypothesis  for  F2  in  the  sense  that  we  now  assume  that  the  time  dependence 
of  F2  resides  within  a  functional  dependence  of  F^ . 

In  line  with  this  hypothesis,  we  expand 

F2(Xj,X2;Fj)  =  F2(Xj,X2;Fp  +  gF2(Xj,  X2;Fj)  +  .  .  . 

and  introduce  operators  that  replace  8Fj/8t  with  in  (96)  and  then  we  equate  equal 
powers  of  g,  for  the  first  two  powers,  to  obtain 


and 


^1^2^*r  *2'^1^  ^0^2^*!' *2'^1^  "  [^2''’'^12'^2^*1' j 


(102) 


■k  ■< 


(2ir)' 


Fj(P^)  L^(ik-Vj) 
1  +  L^(ik-Vj) 


F“(Xj.X2)  L3.(ik-Vj) 


1  +  L^(ik-Vj) 


+  e 


ik-(Q2-Qi)  /  ^  . 

r  J  VPyPpk)  dPj 
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J-.J.UK'  V,) 


1  +  L^Uk- Vj) 


- IT - ^  -^F2(Xj,X2)] 

8F^(Pp8P2V  81*1  / 


dP. 


L^(lk- Vp 

aF°(Xi,X2) 

1  +  L^(ik- Vj) 

aP2 

+  (1 


2)). 


(103) 


Let  us  first  discuss  the  solution  to  (102).  By  the  same  procedure  that  we  used  in 
section  3.  1  we  can  replace  the  D  operator  and  obtain 


By  (26)  the  solution  to  this  is 

F“(x,.X2;si^Fi).sf,F°U,,X2lF,) 

or,  letting  F^— s|^F^  in  the  functional  dependence,  we  have  the  solution 


(104) 


Since  the  left-hand  side  is  independent  of  t,  we  may  again  pick  the  value  of  x  arbi¬ 
trarily.  We  pick  T  —  00  and  use  the  boundary  condition  that  is  similar  to  (49) 


F°(Xj,X2:Fj) - .►Fj(Xj)  Fj(x2).  (105) 

|qi-q2|-oo 

Using  these,  we  obtain 

F°(x..x  :Fj)  =  lim  S^^sJFj(x^)Fj(x2) 

T-^00 

=  Fj(Xj)  Fj(X2).  (106) 


For  the  homogeneous  plasma 

F°(Xi,X2;Fi)  =  Fi(Pi)Fi(P2). 

By  the  definition  of  the  Dj  operator  we  have 
DjF°=DiFj(Pi)F^(P2) 

+  nFj(Pj)  ^^+23(Q2”**3^'^2^^2’F3’^2”‘*3*j  ***3 
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^  2. 
4ire  1 


I ‘">3 * i‘— ^1- 


Upon  inserting  into  (103),  we  see  that  (107)  just  cancels  the  two  integral  terms 
that  have  no  (Qj-Q2)-dependence.  Tliis  cancellation  is  the  plausible  way  of  avoiding 
the  secondary  divergence,  mentioned  with  (19),  which  results  from  an  integration  over 
q^.  The  subtraction  of  DjF2(Xj,  x^)  involves  the  difference  of  two  integrals  that  are 
to  be  cut  off.  This  removes  the  dependence  on  the  cutoff  distance  in  the  equation  for 


o  z 


Thus  (103)  becomes 


^0^2^*!’  *2’^1^  ~  ^2 ’*’'^12' 


4.e2  .  ^Fi<I^2> 

-  n - T  1 - - 


4.e^  ■  ^^l<Pl) 

(2Tr)^  ^  8Pj 


^  ik.(Q2-Qj)  p 


y  fz^Pj.Fl.k)  dpj 


dk.**  ik‘(Q^-Q2)  p  1  ^  ... 


y  ^2<P3'^2'^)  *^3 


p2  ^  p2 

=  Fj(XpX2;F^) 

+  n  y  [4<23(Q^3):Fi(P2)4(Pi.  P3.  07^3)]  (108) 

where,  in  the  latter  form,  we  have  written  the  integrals  in  q-space  and  used  the  fact 
that 

o  +  1^2 

**2  '^IZ  ^  2m  • 

Now  add  the  term  subtract  it  from  the  right-hand 

side  of  (108).  By  comparing  the  results  with  (97)  and  using  the  invariance  of  the  Poisson 
brackets  with  respect  to  a  contact  transformation,  we  see  that  F^  can  be  split  into  two 
functions 

^2(*1'*2’^1^  ^  4(^1’ ^2’^r®2^  ^2^*r*2’^l^’ 

where  h2  is  the  solution  of 

^0*^2^*!'  *2'^1^  “  [ll2^')'l2’^2(*l*  *2’^l^j  ”  ['('l2(^l”^2^’^l(^l^  ^l(^2^j'  (HO) 
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By  the  familiar  use  of  the  operators  we  can  see  that  the  solution  to  (110)  is 

~  00  ^  ~ 

h2{Xi,X2;Fi)  =  -[3^  dT  4.i2(QrQ2-<VrV2)T):Fj(Pj)  F,(P2)J.  (HI) 

Combining  (106),  (109),  and  (111),  we  obtain  a  final  form  for  F^  to  first  order  in  g  imder 
the  adiabatic  hypothesis. 

F2(Xi,X2:Fi)=Fj(P-^)Fi(P2) 

+  fl(Pj,  P2.  ^2)  -  [SJ  4^12<^2-<'^2)^)'Fj(P,)  Fj(P2)  . 

(112) 

In  Section  IV  we  shall  see  that  use  of  this  function  results  in  no  divergence  in  the  equation 
for  Fj. 

We  shall  look  at  (101)  in  its  more  general  form,  as  a  kinetic  equation  for  F2.  As 
such  it  is  exact  in  4*22  accurate  to  first  order  in  all  other  interactions.  If  any  col¬ 
lision  term  in  the  equation  for  F^  is  going  to  be  used,  it  should  be  obtained  from  (101), 
For  example,  the  expansion  in  g  and  the  adiabatic  hypothesis  led  to  (112)  and  will  be 
used  in  such  a  collision  term. 

We  can  also  use  (101)  to  study  the  limitations  of  making  the  adiabatic  hypothesis  for  , 
If  in  the  solution  we  had  expanded  in  g  but  not  made  the  hypothesis,  we  would  have  ' 
had  as  a  zero-order  equation 

-F^^l 

at  ■  [“2^^12’^2j' 

12 

This  is  exactly  the  form  of  equation  that  Bogoliubov  considers  in  showing  that,  for 
short-range  forces,  the  adiabatic  hypothesis  is  valid.  He  shows  that,  for  forces  of  reinge 
r^,  F2(Xj,  X2;t)  would  relax  to  the  form  F2(Xj,  X2:Fj)  in  a  time  of  the  order  r^/v.  If 
for  the  long-range  interaction  we  interpret  a  collision  as  described  in  section  1.  2,  this 
would  be  a  time  of  order 


where  ^  is  the  characteristic  time  of  change  of  Fj^.  We  can  restate  this  by  saying 
that  an  adiabatic  hypothesis  for  F^  is  forced  on  us  by  the  boundary  conditions.  To  zero 
order  in  g,  the  only  force  producing  correlated  motion  is  the  two-body  force,  and  the 


time  that  is  such  that  this  force  is  not  of  order  g  is  small  compared  with  the  time  of 
change  of  F^.  (These  statements  for  F^  can  be  carried  over  immediately  to  the  inho¬ 
mogeneous  case;  but  we  shall  not  carry  this  out  since  the  transition  is  not  as  simple  for 

fJ.) 

If  we  use  the  value  of  F^  thus  obtained  in  the  equation  to  first  order,  v/e  have 
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(113) 


3F2(x  ,X2)  r  ,  -I 

+  "  J  ['*'l3<Qr^3>'^l<^l)  f2<^2'P3’Q?^3>]  ^^3 

+  n  y  [>|-23(^3):Fi(P2)  ^(^l’ P3’ Qr^3>]  "^3- 

Again  an  adiabatic  hypothesis  is  attractive,  since  the  whole  source  term  is  a  functional 
of  Fj.  If  we  call  this  source  term  <j>  (x^,  X2;t),  we  can  write  the  solution  to  (113)  as 

F2(Xj,X2;t^) 


+  ^  S^,_^  <|>(x^.X2;t')  dt' - ►j'  dt'  (|>(Xj,X2;t-t'). 

Here,  we  use  the  boundary  condition  that  F^  is  zero  at  infinite  separation,  goes  to 
zero  for  |q^~<l2l  ^  ^d>  t* -integration  will  back  the  particles  off  to  this  distance 

in  a  time  ~  Xq/v,  If  ^  has  not  changed  in  this  time,  we  hold  it  constant  throughout  and 
say  that 


^00  - 

F2(Xj,X2;t)  =  j  dt'  <t>(Xj,X2;t). 


(114) 


Since  t  in  occurs  only  inside  F^,  (114)  is  equivalent  to  the  adiabatic  hypothesis.  The 

argument  connected  with  (40)  indicates  that  for  a  uniform  plasma  F^  changes  sufficiently 

slowly  that  these  approximations  in  the  integral  of  (114)  are  valid. 

From  this  discussion  we  see  that  the  adiabatic  hypothesis  depends  on  the  order  of 

g.  To  zero  order  it  is  well  founded.  To  first  order  two  problems  enter.  The  first  is 

the  necessity  of  assuming  a  slow  variation  of  F,,  as  mentioned  above.  The  second  is  that 
1 

the  hypothesis  for  F2(x^,  X2)  depends  upon  the  hypothesis  for  F2(x^,  x^)  where  i  #2.  (The 
fact  that  the  source  term  in  (113)  is  a  functional  of  F,  depends  on  this.)  Dupree^ ^  has 
shown  that  the  hypothesis  for  F2(Xj,Xj)  depends  also  on  the  slow  variation  of  F^.  It 
simply  shows  that  the  hypotheses  for  the  various  F^  are  interrelated  and  cannot  be  con¬ 
sidered  separately.  The  first-order  effects  in  a  plasma  include  shielding,  which  is  a 
cooperative  phenomenon . 
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IV.  DISCUSSION  OF  RESULTS 


In  this  section  we  discuss  the  result  obtained  in  the  previous  section  and  see  that 
it  removes  the  divergence  in  (17),  as  desired.  Section  4.  1  will  be  devoted  to  a  dis¬ 
cussion  of  the  equilibrium  case  as  an  example  to  increase  the  understanding  of  the 
results.  Section  4.  Z  will  contain  a  discussion  of  the  more  general  case. 


4.  1  THE  EQUILIBRIUM  CASE 

At  equilibrium  the  momentum  distributions  are  Maxwellian.  Using  the  fact  that 
p  and  P  are  related  by  the  conservation  of  energy, 

-,2  ^Z  2  2 

Pi  *  P;  , 

15^1 ' 


we  may  write  the  first  term  of  (112)  as 


F^(P 


pf  (P^)  =  FfiPy)  Pf (Pz)  exp 


kT|q^-q, 


(115) 


where  F ,  (p)  is  the  Maxwellian  distribution. 

*  Q  10  12 

Other  authors  ’  have  shown  that  for  equilibrium  the  large- separation  solution 
reduces  to 


fJ(Pi.P2-qrq2)  = - ® 


kTjq^-qj 


exp 


Iqrqzl 

Xt. 


Ff  (Pj)  Ff{p2). 


(116) 


In  terms  of  the  undeflected  variables,  and  using  (115),  we  obtain 


fJ(Pl,P2.Q,-Q,)  =  - 


'1 


kT|Qj-Q2 


exp 


iQi-Sj 


Ff  (Pj)  expi 


kT|qj-q2l 

(117) 


For  equilibrium  the  final  term  in  (112)  may  be  evaluated  as  follows: 
[J"  dr  f“(Pj)] 

■  ■  f“(P,)  f“(P2) 

‘  wl”  ^+l2<v^^^|’>  <'Pi'  '■rn’21 
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Ff  (p  )  F^(P2)  exp[ - 

I-Q2I  ^  ^  ^  ^  L 


kTl^l  LkT|q-^|J 

The  results  (115),  (117),  and  (118)  give,  for  the  equilibrium  function  F2(x^,  X2) 
.  „M,_  ,  „M,_  .  _.._r  i  f,.  f,  _ \ 


F^(x^,X2)  =  F"^Pj)  F"‘(P2)  exp 


kT|qj-q2 


kT|Qj-Q2 


1  -  exp 


Those  regions  for  which  IQJ-Q2I  differs  from  |q^-q^  |  are  those  for  which  |q^-Q2|  «  X.^ 
and  for  which,  therefore,  the  last  term  in  parentheses  is  small.  To  the  same  order  of 
eVkTXjj  we  may  replace  |Q^-Q2j  with  jq^-q"^  |  in  (119)-  This  substitution  is  consis¬ 
tent  with  our  assumptions  concerning  the  insensitivity  of  F^  and  F^  to  the  difference 
between  Q  and 

We  could  approach  the  equilibrium  case  from  a  straightforward  application  of  sta¬ 
tistical  mechanics.  It  has  been  shown^  that  the  effective  potential  field  of  a  charge 
imbedded  in  a  fluid  of  charged  particles  is  approximately  the  Debye  potential 

2 


If  we  were  to  assume  that  this  is  the  potential  energy  between  two  charges,  we  would 
write  the  two-body  distribution  function  as 


2 

e  exp 


f|(x^,X2)  =  Ff(p^)  exp  -- 


kT  qj-q2l 


If  ki-qzl  »  e  /kT,  we  can  expand  the  exponential  and  obtain  the  result  corresponding 
to  the  large-separation  solution,  (II6), 


f|(x^,X2)  =  Ff  (pp  Ff  (P2)^  1  - 


kT|qj-q2l 


exp - r - 


If  we  do  not  wish  to  assume  that  »  e^/kT,  we  can  write  as 


F2(Xj,X2)  =  F^(Pj)  F^(P2)  exp  - 


and  expand  the  exponential  of 


kTkj-qj 


kT|qi-q2| 


ki-qj 


kTki-qzl 


r  IV'IzI 

exp - r - 


-  1  . 
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This  function  has  its  maximum  at  jq^-q^l  =  0,  the  value  at  which  it  equals  e^/kT\jj, 
our  small  expansion  parameter.  Under  this  expansion  we  have 


2 

e 

f 

•  4 

kT|qj-4^| 

kT|q^-q2l 


exp 


-  1 


(121) 


Equation  121  is  (119)  when  we  replace  IQJ-Q2I  with  | ,  and  is  the  equilibrium  func¬ 

tion  that  we  obtained  by  allowing  for  possible  small  values  of  jq^-q^j .  Both  (119)  and 
(121)  result  from  expansions  in  the  weakness  of  the  forces  giving  rise  to  shielding. 

If  we  regroup  (121)  as 

r 


F^=Ff(p^)Ff(P2)^ 


exp 


kTki-qzl 


1  +- 


kT|qj-q2l/  kT|qj-q2l 


exp 


and  expand  for  »  e^/kT,  we  obtain  (120)  to  first  order  in  e^/kT|q^^^^| .  The 

effect  of  the  additional  function  h2(x^,X2)  of  section  3.  2  in  the  transition  to  large  sepa¬ 
ration  is  to  cancel  the  first-order  effect  in  the  expansion  of  F^*  We  shall  see  this  in 
the  more  general  case  in  section  4.  2. 

Another  use  for  F2(Xj,  X2)  will  be  in  the  equation 


9Fj(x^) 

at 


Hi;F^(xp 


4^j2:F2(Xi.X2) 


dxo 


(122) 


The  use  of  F2(Xj,  X2)  in  the  interaction  term  of  (122)  involves  the  removal  of  the  dis¬ 
tinction  of  particles  1  and  2,  which  was  discussed  at  the  beginning  of  section  2.  3.  The 
coordinate  X2  in  (122)  is  simply  a  dummy  variable  in  the  integration  over  a  function  that 
is  considered  correct  for  all  values  of  (qp^). 

We  can  use  the  expression  (121)  for  F^  in  order  to  estimate  the  contribution  of 
various  terms  in  F2  to  the  integral  in  (122).  We  know  that  the  angular  integrals  over 
F2  will  give  zero  contribution,  but  we  are  justified  in  examining  the  -depen¬ 
dence  for  this  case.  In  order  to  have  convergence  in  the  -integration,  it  is 

necessary  to  subtract  the  effect  of  the  uniform  background  (also  zero  on  the  angular 
integration).  Thus  we  shall  examine  the  -dependence  of  the  integral 


12 


;F®(x 


r  *2 


)-Ff  (Pi)  Ff  (P2) 


dx2. 


(123) 


Our  motivation  for  this  is  as  follows.  In  the  nonequilibrium  case  the  angular  inte¬ 
grations  of  (123)  will  not  give  zero.  However,  we  can  expect  the  -  dependence 

to  be  roughly  the  same.  If  we  can  solve  the  equilibrium  case,  we  can  gain  some  in¬ 
sight  into  the  nonequilibrium  case. 

Using  (121)  for  f|  in  (123),  we  can  write  the  part  of  the  dx2- integration  as 
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(124) 


P  00 

r  ^ 

/  2  \ 

2 

r  2 

\  dr 

*^0 

1 

1 

1 

_e _ r_ 

kTr  Xjj 

(124) 

15 

Let  us  integrate  the  last  term. 


2  n  ( 

w  L 


dr 


exp 


kTr 


We  can  expand  the  modified  Bessel  fxmction  for  small  argument  and  obtain 


2 

e 


kTr 


(125) 


C 

where  y  =  e  ,  and  C  is  Euler's  nximber,  0.  577. 

Use  of  the  function  (121)  has  resulted  in  a  convergent  integral.  Had  we  followed 

other  authors  and  used  (116).  the  integi‘'al  would  have  diverged  at  small  distances. 

2 

Cutting  it  off  at  the  distance  of  closest  approach,  e  /kT,  yields 


2  p  00 

-^e^/kT 


2 

This  result  agrees  with  the  more  exact  treatment  to  the  same  order  of  e  /kTkj^. 

This  surprising  accuracy  of  the  cutoff  procedure  occurs  because,  to  a  function  that 
varies  as  the  function  exp(-e  /kTr)  cuts  off  very  sharply.  Because  this  result  was 
obtained  only  from  the  -dependence  of  the  integral  of  (123)  and  did  not  rely  on 

the  symmetry  of  the  equilibrium  function,  we  can  expect  similar  results  for  any  non- 
equilibrium  function  f2(P2,  We  therefore  expect  that  we  can  still  retain 

accuracy  to  the  same  order  in  g  if  we  set 


4(Pi.P2,Qi-Q2)  = 


*^2^^1  ’  ^2'  ^1  ~^2^ 


(126) 


2 

and  cut  off  the  integration  at  e  /kT. 

The  new  results  obtained  by  considerations  of  the  close  approach  are  in  the  first 
two  terms  of  the  integral  in  (124), 


(127) 


This  is  of  the  same  order  of  g  as  (125),  and  thus  we  obtain  a  non- negligible  contribution 
as  a  result  of  the  close  collisions.  The  generalization  of  (127)  to  the  nonequilibrium  case 
will  not  be  as  simple  as  the  result  for  f^,  (126).  The  generalization  of  (127)  depends 
very  intimately  on  the  nature  of  the  functions  for  |q2-q^  |  ~  e^/kT,  and  will  be  inves¬ 
tigated  in  the  next  section. 
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4.2  THE  NONEQUILIBRIUM  CASE 

The  discussion  relating  to  the  general  case  is  necessarily  more  complex,  since  we 
must  deal  with  general  expressions  rather  than  evaluate  integrals  of  known  functions. 
We  shall  see  that  to  first  order  in  g,  the  entire  effect  of  the  close  collisions  can  be  put 
into  one  term  that  is  a  generalization  of  the  Boltzmann  collision  integral.  To  this  will 
be  added  another  accounting  for  the  velocity  dependence  of  the  collective  interaction. 

Consider  (99)  for 


fJ(Pl.P2.k)  = 


4,rie2  k  Di2Fi(Pi)  F^(P2) 


-  k 


ik-(Vj-V2))L  aPj  ^  ^  ^  ^ 

<$3  f‘(P,.P3,k)|. 


In  (128)  we  define 

f5(Pj.  P2.  k)  =  ^  d^Q  e^  ■  f^iPj.  P2.  Q) 

and 


(128) 


(129) 


k^  (€-it-(7;^)) 


'*'12<Qi2-<Vi-V2)t);Fj(Pj)Fj 


■l(?2)]. 


(130) 


From  (130)  we  see  that  smeill  distances,  IQJ-Q2I  «  Vj^,  correspond  to  large  wave 
numbers,  k  »  in  the  region  of  chief  contribution  of  the  transform.  For  l2u:ge 

k,  (129)  can  be  expanded  in  inverse  powers  of  k  by  using  the  first  term  of  (129)  as  the 
first  approximation  in  an  iteration  expansion. 


nf^(Pj.P2.k)+^i 


(€-ik(Vj-V2)) 


(€-ik-(Vj-V2)) 


r  ^•P32^1<P3) 

8Pj  k.(V3-V2) 
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-  8F  (P,) 

-  k - 

aPo 


P2  y  k-<Vj-V3)  J 


Successive  terms  fall  off  as  k 

For  k  >  ykr/e^XQ  the  right-hand  side  of  (131)  is  of  order  {e^/kTX.jj)^,  which 
is  of  second-order  smallness  in  g.  As  was  pointed  out  in  section  1.2,  the  differ¬ 
ence  between  P  and  p  and  Q  and  q  becomes  important  only  in  the  region  |q,-qol  < 
r  •>  A  ^  ^ 


ye  \j^/kT.  This  fact  implies  that  on  the  left-hand  side  of  (131)  we  can  replace  Q 
with  q  and  P  with  p  and  know  that  in  those  regions  in  which  this  replacement  is  not 
valid  the  function  (131)  is  very  small.  Thus  in  q-space 

fJ(Pl.  P^.  ^2)  -  [ ^  “  dT  ^j2(QrQ2-(^^2)^)'Fj(Pp  F^{%)\ 

- *"f2(Pl-P2- dT,j-(qp:t^-(7^2)'"):Fj(p;)  Fj(r2)].  (132; 

where  the  error  incurred  will  be  of  order  (eVkTXp.)^.  (The  case  of  |q,  -Q-  |  ~ 

V2  ~  ~  U  1  c 

e  Xjj/kT  will  be  mentioned  later.) 

The  final  use  for  F2(Xj,  X2)  will  be  in  the  interaction  term  (I.T.)  of  the  equation 
for  F^, 

»Fi(x  )  r  -1 

— i^-[H°:Fj(xpJ=I.T. 

»  n  y  [^^12’'|f^(P^)  Fj(P2)+f2(Pi>  P2> 

-[Iq  dT  ^>j2<^r^-<^^>^)’Fj(ri)  F^(^)]}]  dx2.  (I33; 

where  we  have  proved  that  there  is  no  divergence  for  |q^-q^ |  —  0. 

The  third  term  of  (133)  appears  to  diverge  for  |qj-q^  |  —  00,  but  we  can  easily 
show  that  this  divergence  is  cancelled  by  the  first  term.  Consider  the  equation  for 
Fj(Pj)  Fj(P2)  in  the  region  of  small  i|<j2>  that  is,  large  , 

^0^2^*!’  *2^  ~  [^2’F2^^1’  *2^]  [^12’F2^*1’  ^^2^  ‘ 

Using  the  usual  solution  methods  and  (26),  we  may  obtain  an  expansion  for  small  t|«^2> 
Fj(Pj)  Fj(P2)  =  Fj(?^)  Fj(r2) 


^  [L  +12<‘>l-^2-<^l-'^2>^)  dr;  Fj(ri)  Fj(?^)] 
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+  ...  .  (134) 

Using  this  expansion  in  (133)  in  the  region  of  large  |qj-^| ,  we  see  that  the  second  term 
of  (134)  cancels  the  divergent  term  of  (133).  The  higher  powers  of  i|'22  (134)  will  con¬ 

verge  in  the  integral  of  (133)  for  (q^-q^l  —  oo. 

We  see  that  (133),  although  a  sufficient  result,  is  in  an  inconvenient  form  because 
the  final  term  in  the  interaction  term  cancels  the  divergence  in  fg  as  |qp^|  —  0  and 
the  divergence  in  Fg  as  Iqj-q^j  -►  oo.  We  shall  obtain  another  form  that  is  equivalent 

to  (133)  but  that  is  more  convenient  in  the  kinetic  equation  for  F,. 

^  1 

In  the  Appendix  we  derive  the  result  (17)  corresponding  to  the  f2  term  in  the  inter¬ 
action  term  of  (133), 


dr 


I6v2ne^ 

(2.)’  JJ  k" 


a  6(k-(Vj-v)) 

- - 5  k  •  d. 

ap  |l+L  (ik-v  )p 


Fi(Pi)  Fi(P2)' 


(135) 


The  corresponding  expression  for  the  last  term  in  (133)  is  obtained  by  expressing  the 
integral  in  k -space  and  is  given  by 


,/  2  4  dkdpp  r,  _  ^ 

"  ”3  \\— ^k--^6(k-(v^.))k-d  F  (p  )F  (?).  (136) 

{Zirr  k^  apj  ^  ^  1  l  1  ^ 

(Using  the  fact  that  L^(ik  •  v)  0,  we  can  easily  see  that  the  difference  between 

(135)  and  (136)  converges  for  large  |k|  in  agreement  with  (132).) 

In  the  interaction  term  of  (133)  add  and  subtract  the  term 


I6v!neirr_^^ 

(2ir)^ 


6(k.(vp=r2))k-dj2Fj(Pj)Fj(r2) 

api 


dT  >|>j2(qi-q2-(Vi-V2)T):  F 


l(p,)Fi(r2)] 


dx, 


2’ 


(137) 


where 


and 


iD  _ 
'*'12  = 


hi-q 


—  exp 


qrq. 
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The  difference  between  (135)  and  (137)  also  converges  for  |k|  -  oo. 

With  these  additions,  the  interaction  term  becomes 

I.  T.  =  n  y  [+^2=  F^(Pj)  Fj(ri)  F^(P2)]]  dX2 

^  "  I  Fl<Pi)  Fi(P^)]]  dX2 

J{[‘*'l2’  ^Z^Pl’Pz’ ‘^1~‘^2^] 

■  Fi(Pi)]]}  dx2.  (138) 

In  order  to  combine  the  first  three  terms  of  (138),  consider  the  equation  for 
^l(Pl)  Fl(P?)  '-here 


P? 


j  =  lim  st;pj. 


P?  =  lim  SD  P2. 

T-*00 


(139) 


and 


_a_ 

3t 


for  any  <Kx^,  X2).  Here,  is  the  operator  that  projects  the  particles  1  and  2  back 
along  trajectories  corresponding  to  a  Debye -shielded  collision. 

By  the  familiar  means,  we  have 


Vl(f?)  =  [h?<^F,(p5>)  F,(pf)] 

'  [h?**.,.f.(p?)  f.(f?)] 

*[*?2-l'irF,(p?)  p,(pf)]. 


(140) 


Since 


I'*'l2  '^izl  /^D’ 

we  can  expand  Fj^P^^  F^^P^^,  using  (26),  in  terms  corresponding  to  powers  of 

(eVkTXjj), 

'■l(^)  -  F,(P,)  F,(P2)  .  dT  -|.,2iF,(P,)  F,(P^)].  ,141) 
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The  difference  i)>^2  ~  ^12  very  small  imtil  the  r-integration  backs  the  particles  to  a 
separation,  ~  X.^,  so  that  we  obtain  in  the  integral  of  (141)  only  the  asymptotic  velocities, 


(QjH52-(V,-V2|r)-+,2lQi-Q2-(Vi-Vj)T)};  F,!?,)  Fj(P2l].  (142) 

As  in  (132),  this  difference  is  small  unless  jq^-qg  |  5  kjj,  in  which  case  Q  -►  ^and  P  —  p 
(|Qj~Q2I  ~  ^£)  discussed  later^.  Thus  (141)  becomes 

".(*?)  *  ''i(^)  ''1(^2) 


(143) 


Using  (143),  we  can  obtain 

I  hz:  Fi(P;)  F^(P;)]  dX2  2  J  {  [4<?2^  F^(P?)  Fj(p?)] 

-[‘*‘?2'  [Iq  +?2<^-<^2)^)'  Fl<^l)  ^1<^2)]]}  ^2' 

(144) 

2 

where  the  equality  holds  to  first  order  in  (e  /kTXp).  Using  (144)  in  (138),  we  have 
I.T,  =  ny[il.°:F,(p°)  F^P?)]  dx2 

J  {['•'iZ’^Z^Pl'Pz’ ‘ll"‘*2^] 


We  may  see  from  the  following  considerations  that  the  first  term  of  (145)  is  a 

Boltzmaiui  collision  integral  with  a  Debye -shielded  interaction.  The  method  used  is  due 
12 

to  Bogoliubov. 

From  the  definition  of  p° 

[hJ++“  lF,(gf)  F,(^)]  •  0. 

Thus  we  set 


I  [+?2‘'’l(^)''l(^)]  <^2  '  ''l(^) 
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*D 


and  recognize  that  depends  on  only  through  q  =  q2  -  For  the  q- integration 
set  up  a  cylindrical  coordinate  system  with  the  positive  z-axis  parallel  to  V2  - 
Denote  the  radius  and  polar  angle  by  a  and  4>  >  respectively.  Then 


^  c  r ""  4! « A w) 


(147) 


-^D  '^D 

From  the  two-body  problem  under  consideration,  (Xj,  X2)  and  P2  (x^,  X2)  are  the  pre- 
collision  momenta  that  will  yield  the  state  Xj,  X2.  From  this  definition 


Pi  (Xj,  X2) 


P°(Xi,X2) 


4=-00 


4=+00 


=  p 


1 


Pi 


P2  (Xi,X2) 


P2  (Xi*  X2) 


i=-CO 


^=+00 


=  Po 


(148) 


Pz 


where  Pi  and  P2  are  those  momenta  that,  as  precollision  momenta,  will  give  Pi  and  P2 
as  a  final  momentum.  For  coordinates  x,  and  x.  which  are  such  that  an  interaction  is 
in  progress,  P.  and  p.  are  not  equal.  Only  when  the  particles  are  separated  by  a 
distance  that  is  much  greater  than  do  PV  and  Pi  become  equal. 

Using  (146)-(148)  in  (145),  we  have 


♦lo  ‘'0 


P7-P 


(P2> 


(149) 


which  is  a  form  of  the  Boltzmann  collision  integral  with  the  Debye-shielded  interaction. 
Using  (135),  (137),  (145),  and  (149)  in  (133),  we  obtain 


®^l<*l’  r  o  1 

— ir^=[H°;Fi(Xi)J 


IPz-pJ  (Fl(^)^l(P2)-Fl(^)Fi(Pi))  dp^adad* 


(P2> 


.  I6ir^ne^  f  r- 

+  - 5 —  \  dk  dp,  k 

(2it)-*  ^ 


._l/ _ L_ _ 

8p^\k^|l+Li.(ik  - Vj) 
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(150) 


'  o?Tir7j  ^ 

6  ••  2 

The  collision  integral  is  the  same  as  that  used  by  Allis  and  by  Spitzer  and  Harm 
insofar  as  they  replace  the  Coulomb  field  with  a  Debye  field  in  adapting  the  Boltzmann 
theory  to  plasmas.  It  is  convenient  that  the  entire  effect  of  close  collisions  is  contained 
in  this  term.  This  is  consistent  with  the  concept  of  the  close  collision  and  its  analogy 
with  the  Boltzmann  gas  discussed  in  section  1.1. 

The  integral  over  k  in  (150)  is  the  contribution  to  the  interaction  term  which  is  due 
to  the  deviation  of  the  shielding  cloud  from  the  Debye  shield.  This  is  caused  by  the 
nature  of  Fj  and  the  velocity  v^,  as  can  be  seen  from  the  fact  that 


lim  L  (ik  •  ^, )  =  lim 

Vi-0  ^  v^-0 


2 

moj 

P 


I 


k  •  (vp^) 


9F^(u) 

9u 


du 


u 


(151) 


B’or  a  distribution  that  is  spherically  symmetric  the  integral  is  real.  For  a  Maxwellian 
distribution  (151)  reduces  to  k^/k  ,  which  makes  the  k  integral  of  (150)  zero. 

If  the  distribution  Fj(p)  is  not  Maxwellian,  the  contribution  of  the  integral  over  k 
in  (150)  is  of  the  Fokker-Planck  form.  If  we  define  D°  and  B°  in  (5)  we  have 


jjO  =  _  I6ir  ne 
m(2Tr)^ 


dp2dkk  j 

k" 

ll+L^dk-iTj)!^ 

V 

1 


.  *  _ _  ^  9F  (p^) 

■r  6(k.(v^))  k - 

zrl  ^  ^  apl 


1  +- 


rr  dr^dkkk 

"  2,,_.3  JJ  .4 


k 

.  \ 


(152) 


m^(2ir)'*  k^  I  |l+L^(ik- Vj)|  ^ 


rfjj 


6(k  -  (v^-Vj))  Fj(P2). 


Armed  with  solution  (150),  we  are  in  a  position  to  justify  the  assumptions  (90), 
(95),  (132),  and  (142).  We  have  shown  explicitly  that  these  assumptions  are  good  for 

IQ1-Q2I  <<  ^ 

We  can  now  show  that  the  error  introduced  by  these  assumptions  for  Iq^-Q^I  ~  X.^ 
is  negligible.  We  are  interested  in  an  integral  of  F2(Xj,X2)  over  and  In  this 
integration  consider  a  value  of  (q^-tf^)  which  is  such  that  Iq^-Q^  |  ~  (which  also 

means  that  Iq^-^l  ~  X^).  The  number  of  particles  at  this  separation  which  have  under¬ 
gone  a  close  collision  in  the  past  is  given  by  the  number  of  particles  that  would  undergo 
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a  close  collision  in  an  inverse  collision  starting  from  the  same  spacial  coordinates  with 
the  velocities  reversed.  The  probability  for  this  is  given  by  the  solid  angle  that  is 
available  for  close  collisions  and  that  is  given  by  A/R  ,  where  A  is  the  cross  section 
for  close  collision  and  R  is  the  initial  separation  of  the  particles.  For  particles  sepa- 
rated  by  this  ratio  is  (e  \j^/kT)/\j^  =  e  /kTX.Q.  Since  nothing  radical  happens  for 
these  particles,  for  example,  we  do  not  predict  a  divergent  result  for  their  contribution, 
we  may  say  that  the  approximations  (90),  (95),  (132),  and  (142)  are  valid  for  all  those 
situations  that  yield  a  significant  contribution  to  the  kinetic  equation  for  F^. 

As  a  final  step  we  shall  evaluate  the  )kj  -integration  for  B®  and  D®  in  (152),  which 
will  leave  these  functions  expressed  in  angular  integrals.  The  |k|  -integration  can  be 
performed  in  the  general  case  as  follows.  Define 


L^(ik-^j)  =  k^L^(ik-T^) 


(153) 


and  note  that  does  not  depend  on  |k| ,  only  on  the  direction  of  k/k. 


W(0,<(.,Vj)  = 


We  define 


(154) 


where  0  and  ^  are  the  polar  and  azimuthal  angles  of  k  in  k- space.  W  converges  at 
the  upper  limit  as  a  result  of  the  subtraction  of  the  two  terms,  and  the  integral  in  (154) 
can  be  evaluated  by  standard  means  to  yield 


W(0,<|).vj)  =  Re 


l;  -  l;  log*  -f 
_ ^ 


l;-li 


(155) 


In  terms  of  W,  the  functions  B°  and  D°  can  be  written 


— “  1  / _2  4  p  p.2ir  P If  ^  ®F,(^) 

dO  =  -  \  dm  \  d«(>  \  sin0  d0  ■^W(0,.).,r,)  6(k-(v^))  k - 

m(2Tr)^  ^  -^0  -Jq  ^  ^  ^  8p; 


m(2ir) 
2  4 


(156) 


For  a  spherically  symmetric  distribution  Fj(p),  it  is  convenient  to  pick  the  z-axis  of 
the  k- integration  along  the  direction  of  For  this  case  the  azimuthal  integration  yields 
2ir,  and  B°  and  are  expressed  as  one- dimensional  integrals.  The  final  evaluation 
of  these  integrals  must,  of  course,  be  made  with  a  particular  F^(p). 

In  the  use  of  (150)-(156)  the  following  should  be  noted.  Although and  are  in 
the  form  of  ^''Tdter- Planck  coefficients,  they  do  not  represent  the  complete  contribution 
to  the  dynai  il  friction  and  the  diffusion  in  velocity  space.  One  method  of  treating 
the  Boltzmann  collision  integral  is  to  expand  the  integral  in  a  series  of  powers  of 
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2  f> 

deflections  in  velocity.  ’  The  coefficients  of  the  first  two  powers  of  this  expansion 
will  give  additional  contribution  to  the  dynamical  friction  and  the  diffusion  in  velocity 
space.  The  contributions  (136)  are  due  to  the  particular  shape  of  the  shielding  cloud 
and  happen  to  be  in  the  form  of  Fokker-Planck  coefficients. 
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V.  SUMMARY 


Equation  (150)  is  the  result  that  was  sought.  As  a  kinetic  equation  for  F,{p),  it  may 
be  considered  accurate  to  first  order  in  (e  /kT\j^).  The  use  of  the  adiabatic  hypothesis 
has  yielded  the  most  general  kinetic  equation  for  Fjj  and  insofar  as  the  hypothesis  holds, 
this  equation  can  be  said  to  exist. 

The  form  (150)  is  appealing  because,  although  it  is  a  new  result,  it  contains  terms 
corresponding  to  those  of  various  earlier  treatments,  the  Boltzmann  equation  (3)  and 
the  Fokker -Planck  equation  (5).  (There  is  no  Vlasov  type  of  contribution  for  the  uni¬ 
form  plasma.)  Equation  3  was  correct  for  close  collisions  and  in  error  for  grazing 
collisions;  and  (5),  the  opposite.  Equation  150  combines  these  two  solutions  and  forms 
a  bridge  between  them. 

Discussions  of  the  solution  of  (150)  will  not  be  carried  out  here,  since  the  various 

2  ^ 

terms  in  (150)  correspond  to  earlier  treatments.  Spitzer  and  Harm  and  Allis  have 
discussed  the  Boltzmann  collision  integral  with  an  assumed  Debye -shielded  interaction. 
Several  author s®~^^  have  obtained  the  term  containing  the  integral  over  f^  in  (150). 
This  result  was  (17)  mentioned  in  section  1.2.  The  divergence  in  the  k -integration  is 
cancelled  by  the  last  term  of  (150),  although  the  other  properties  remain  the  same. 

The  solution  of  (150)  is  still  complicated  and  must  be  subjected  to  other  approxima¬ 
tions.  However,  within  the  limitations  mentioned  above,  (150)  may  be  taken  to  be  the 
equation  of  evolution  of  Fj^,  and  these  approximations  of  solution  may  be  used  and 
analyzed  in  their  own  right.  The  ability  to  separate  the  approximations  used  in  deriving 
the  equation  from  the  approximations  used  in  its  solution  aids  greatly  in  understanding 
the  accuracy,  meaning,  and  limitations  of  the  analysis. 
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APPENDIX 


We  shall  carry  out  a  derivation  of  the  large -separation  solution  £2  which  is  due  to 
Dupree.  We  shall  then  show  that  Dupree's  result  can  be  put  into  the  form  obtained 
by  Lenar  d. 

•'  ,  1  9  . 

Dupree  considers  (97)  for  q  with  the  substitutions  —S_^  and  D^  — as  in 

(66).  We  write  ^  -  q2  and  obtain 


at  4  ^Pl’  P2’  S_^Fj) 


-n  j  dr3dq'>)-(q-q')f^  (P3.  P^.  q’:  S^^F^) 

9Fi(^)  a  r  ^  1 _  1 

+  - if  3  dp3dq'4-(q-q')f2  (Pi.P3.q';S_^Fj) 


dj2F^(p;)Fj(^). 


(157) 


This  can  be  written 


(158) 


where  Lj  and  L2  are  Landau  operators  defined  in  (71)  and  operate  on  coordinate  1  and 
coordinate  2,  respectively.  Note  that  L^  and  L2  commute. 

The  formal  solution  to  (158)  is 


00  -(l,+L-)t'  a  . 


^v|j(q)  •  dj2  Fj(Pj)  Fj(p2) 


where  we  have  taken  the  limit  t  —  00  in  agreement  with  the  discussion  preceding  (98). 
It  is  easier  to  evaluate 


J  d^  fj(p;.p^,q) 


and  to  then  compute  from  (157).  Using  (78)  and  (79)  for  the  result  of  the 


operator  L^,  we  have 
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\f^  (r..r2.k)d^  = - 4-  \  dT  \\  da  dcr,  firT-TirT  •  \  - 

j  ^  ^  ^  2  Jg  _  J_J  1  2  1  +  L_(  a2)  ^  .JJ-,  ^ 


^  ®F  (^) 

_ P  k - - — — 


_!sf  -, _ !L_  f 

*  T, -ik-r  •>+ 


dp; 


j  .  aj-ik*Vj  aj-ik*Vj 


4irc 


2  '  **12 


(159) 


where 


f2  (P‘i*P^*k)  =  j'dq  e^^‘ ^f2  (fj.P^.  q). 

Equation  159  is  the  form  obtained  by  Dupree.  The  ^  dp  and  are  defined  in  (78) 

and  (79).  The  ^  dp  and  L_  are  the  same  integrals  with  the  path  of  integration  going 
under  the  poles. 

As  mentioned  in  connection  with  (80),  the  a-integration  can  be  moved  just  to  the 
left  of  the  imaginary  axis  and  the  r-integration  carried  out.  This  brings  down  a 
(ai+a2)  ^  if  we  can  assume  Re  a2  <  Re  aj  (if  not,  carry  through  the  following  oper¬ 
ation  with  aj  and  a2  interchanged).  For  the  same  reasons  and  conditions  that  the 
zeros  of  1  +  L^(a)  are  in  the  left-half  plane,  the  zeros  of  1  +  L_(-a)  are  also  in  the 
left -half  plane.  Therefore,  we  close  the  a2  in  the  right -half  plane  and  have  only  the 
contribution  at  a2  =  -a^.  This  yields 

j-r*  dp*  4  d, 


^  ”"«p  ^  8F  (p-j) 

-^ik. 


[l  +L^(a)]  [a-ik-^] 


r  d^^ 

-ik  •  ^  1  *'+  a  -  ik*  ^ 


a  -  ik* ^ 


^  '  ^12  ^1<^1>  ^1<^>  • 

By  algebraic  manipulations  (160)  can  be  put  into  the  form 


(160) 
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(161) 


2  r“i+P  d(r _ 1 _  r 

irk"  ‘^-ooi+p  <r-ik-^  ( 1+L^(<r))  ( 1+L_(a))  J_ 


I  8F  (^)  p  dp; 


In  the  limit  of  very  small  p 

'z  -  - 2 


*<^•^12  ^i(Pi)  ^i(^)- 


‘  f  .,  !L^-.  . J- * ,  f  5- 

“4<r-ik.^  ”^4-ik.^  ^  ^  ™ 


do-  Q  d<r  .  -C.  r. 

\  - =  X - +  wi  \  do-  6(o--ik*v.). 

•^-ooi-^p  o-^ik-v,  ‘^o--ik*v, 


(162) 


Using  (162)  in  (161)  with  further  algebra,  we  obtain 


(  ^2  <^1*  ^*^2  =  I  - - - 

J  ^  ^  ^  2  ^Z  -l-ooi+p  (T -ik-^j  (l-tL^(<r))(l-l-L Jo-)) 


(1+Ljv))  (  - •  ^'^2  ^1^^^ 

*^+  V  -ik.^  14  1  1  14 


2v  J  dp;  6<k-^-t-ia)  k*d^2  ^l^^l^  ^1^^^ 


The  first  integral  may  be  closed  in  the  right-half  plane  with  zero  contribution  because 
there  are  no  poles  in  this  plane.  (By  reflection  the  zeros  of  1  -t-  L  Ja)  are  in  the 
right-half  plane  because  the  zeros  of  1  -I-  LJ-o-)  are  in  the  left-half  plane.) 

We  obtain,  as  a  final  result, 

j* f2*Pi.^.k) 

.±L£!  _ 1 

k^  ‘^-ooi-t-p  <r -ik.^j  (l-»-L^(v))(l-tLJa)) 

Jd^  6(k.^-Wa)  k-dj2  FJ?;)  Fj(5;) 
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r  —^rrrz  *1 - rr  6(k.^+i«r)  k  •  d  F  (fSt)  F  (^) 

<r-ik.^  1+L  1-5  1  1  I  Z 


k  (T-ik*^  1+L^(<r)| 


k'^  I  1+L^(ik-Vj) 


6(k.(v2-Vj))  k.^2  Fi<P2>  * 


where  we  use  (162)  and  the  fact  that  (1+L^((j- ))( 1+L_(ir))  =  |  l+L^(o-)  |  for  o-  that  is  pure 
imaginary. 

Fortunately  the  principal  value  term  in  (163)  is  even  in  k  while  the  second  term  is 
odd.  Thus  in  the  integral 

2. 

mto  1 

j't^l2’  “  “  ^2it)3  .^^2  <Pi'P2'^>  ‘^Pz 

we  need  only  the  second  term. 

Introducing  (163)  into  (164),  we  obtain  (17)  for  the  evolution  of  a  specially  homogene¬ 


ous  Fj, 


®^i<^>  i6Ac^  rr r  &  6(k.(v2-vp) 

~'"5PT  W- 
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